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ABSTRACT 


A simplified  one-dimensional 'analytical  model  for  the 
steady  state  combustion  of  a multiparticle  cloud  of  boron  in 
air  was  established.  This  model  stressed  the  importance  of 
the  radiant  heat  exchange  between  the  luminous  flame  and  the 
burning  particles  adhering  to  the  combustor  walls  with  the 
cold  gas-particle  mixture.  The  model  assumed  no  thermal  or 
dynamic  slippage  between  the  fuel  particles  and  the  carrier 
gas.  Further  assumptions  included  spherical,  uniformly 
distributed,  one-micron  diameter  particles  and  an  optically 
thin  cloud,  during  heat  up  to  ignition  temperature.  The 
particle  reaction  model,  based  on  a single  particle  reacting 
in  air,  was  formulated  to  consider  the  finite  environment  of 
oxygen  and  the  temperature  dependence  of  the  reaction. 

A steady  state,  self -sustained  flame  of  elemental  boron 
in  air  at  atmospheric  pressure  was  investigated  to  evaluate 
the  theoretical  analysis.  A particulate  cloud  of  boron  was 
Ignited  using  a hydrogen  flame;  once  the  flame  was  establish- 
ed, the  hydrogen  supply  was  removed.  The  flame  was  then 
allowed  to  stabilize  as  a self -sustained,  steady  state  flame. 
Combustion  gas  sampling  and  absorption-emission  measurements 
were  taken  to  evaluate  the  combustion  efficiency  and  tempera- 
ture. These  data  were  taken  at  the  combustor  exit  plane  for 
combustors  1.9  to  7.62  centimeters  in  diameter  and  up  to  ?6 
centimeters  in  length.  Particulate  concentrations  providing 
fuel-to-alr  equivalence  ratios  ranging  from  one  to  two  at 
air  mass  fluxes  up  to  6.35  gm/cm2-seo  were  investigated. 


"M)WW> Wfcftw  WWWIWW  . W-»  V.w~. -....  


DS/ME/74-2 

Trona  amorphous  boron  powder,  with  a nominal  particle 
size  of  one  micron,  was  used  as  the  fuel.  A fuel  feed 
system  was  devised  that  permitted  measurement  of  the  boron 
powder  flow  rate  and  provided  a steady  flow  of  fuel 
particles.  A portion  of  the  boron  adhered  to  and  burned  on 
the  combustor  walls  providing  a refractory  surface  that, 
coupled  with  the  energy  radiated  from  the  flame,  produced 
sufficient  energy  to  heat  the  oncoming  particle  cloud  to  . 
ignition  temperature  when  the  chamber  residence  times  were 
in  excess  of  two  milliseconds. 

Complete  oxygen  depletion  and  flame  temperatures  near 
the  adiabatic  flame  temperatures  were  recorded  for  combustor 
lengths  in  excess  of  6l  centimeters  and  combustor  residence 
times  greater  than  five  milliseconds.  Experimental  combus- 
tion characteristics  calculated  from  the.  analysis  yielded  an 
Arrhenius  activation  energy  near  40  Kcal/mole,  an  Arrhenius 
pre-exponential  factor  of  9*04  x 10~3  moles  of  alr/cm-sec 

and  a d-law  evaporation  rate  constant  of  949  /t2/ sec. 

2 

Although  a d-law  or  diffusion  controlled  reaction  analysis 

was  assumed,  these  findings  correlated  best  with  previous 

results  using  a d-law  or  chemical  controlled  reaction. 

2 

However,  the  d-law  analysis  was  used  to  be  consistent  with 


previous  investigators 
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I.  INTRODUCTION 

There  has  been  a continuous  interest  in  the  combustion 
of  metals  for  the  past  few  decades.  This  Interest  has  been 
due  to  the  high  volumetric  heating  values  and  high  combus- 
tion temperatures  possible  with  these  systems.  In  spite  of 
this  interest,  oar  present  understanding  of  metal  combustion 
lags  far  behind  that  of  conventional  gaseous,  liquid  and 
solid  fuels.  The  reasons  for  this  lag  in  understanding  lie 
not  only  in  the  lack  of  research  efforts  on  metals  as  an 

energy  source  but  on  the  very  aspects  of  metal  combustion 

_ * 
that  make  It  attractive  as  a high  energy  fuel.  The  large 

exothermic  heat  of  formation  of  the  combustion  products 
implies  that  they  are  condensed  phase  products  which  are 
highly  stable.  Thus,  the  fuel  and  the  products  of  combus- 
tion are  present  in  the  combustion  zone  in  the  liquid  or 
solid  state  at  normal  combustion  temperatures  (Ref  1), 

Boron,  like  some  other  metals,  has  been  shown  special 
Interest  because  of  its  high  theoretical  energy  content.  Of 
the  basic  elements,  listed  in  Table  I of  Appendix  J,  only 
hydrogen  and  beryllium  have  greater  energy  content  per  unit 
weight  bhan  boron;  but  due  to  their  lesser  densities,  boron 
has  the  highest  heating  value  per  unit  volume.  However, 
boron  as  a propellant  additive  has  been  very  disappointing 
because  it  exhibits  serious  combustion  Inefficiencies. 

Boron  has  been  burned  as  a liquid  slurry  in  a ramjet  (Ref  2) 
and  as  a solid  propellant  additive  in  a ducted  rocket  (Ref  3)» 
The  primary  problems  encountered  in  these  instances  were  the 
difficulty  of  igniting  the  boron  and  the  slow  reaction  rates 
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resulting  in  unreacted  boron  in  the  exhaust  plume.  Many 

theories  have  been  postulated  as  to  the  causes  for  this  poor 

performance.  Primary  among  these  ares  (a)  the  formation  of 

an  inhibiting  oxide  layer  on  the  metal  particles,  (b)  the 

agglomeration  of  tbs  particles,  and  (c)  the  formation  of 

unstable  products  of  combustion  that  inhibit  the  reaction 

and  thus  reduce  the  heat  of  combustion  (Ref  5)* 

Boron  has  other  interesting  characteristics  as  well. 

Technically,  boron  is  not  a metal;  it  is  the  only  element 

not  classed  as.  a metal  that  has  less  than  four  electrons  in 
* * 
its  outer  'shell.  Prom  the  standpoint  of  safety,  boron  . 

differs  from  the  other  metals  in  that  it  presents  no 

problem  from  bulk  fire  hazards,  spontaneous  dust  explosions, 

or  toxic  combustion  products  (Ref  4).  From  thermodynamic 

considerations  at  ambient  conditions,  boron  is  the  only 

metal  whose  adiabatic  flame  temperature  in  air  can  exceed 

the  vaporization  temperature  of  its  oxide.  As  shown  in 

Table  II  of  Appendix  J,  it  is  one  of  the  few  metals  that 

evaporate  at  a temperature  higher  than  that  of  its  primary 

oxide  and  also  well  above  its  adiabatic  flame  temperature, 

approximately  3000°K.  These  characteristics  Indicate  that 

while  BgO^  is  liquid  or  solid,  it  forms  a protective  coating 

on  the  boron  surface.  This  coating  Inhibits  low  temperature 

combustion;  thus,  only  when  the  protective  coating  has 

evaporated  will  the  reaction  rate  of  boron  be  fast  enough  to 

be  a satisfactory  propellant  additive. 
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Background 

The  rates  of  the  ignition  and  combustion  processes  of 
metals  have  been  found  to  be  directly  proportional  to  the 
subdivision  of  the  metal  fuel;  therefore,  most  studies  of 
metal  combustion  have  investigated  fuel  particles  with 
diameters  less  than  100 /u  (Ref  1).  Most  fundamental  studies 
have  been  carried  out  on  single  particles;  although,  the 
qualitative  effects  of  concentration  of  dust  clouds  have 
been  investigated  (Refs  5 through  14).  The  investigations 

with  particulate  boron  have  shed  light  primarily  on  its 

* « 

ignition  characteristics.  The  best  evaluation  of  the 
combustion  characteristics  was  accomplished  by  Talley 
(Ref  15)  using  bulk  boron. 

Talley  studied  the  combustion  of  boron  in  oxygen  using 
a boron  rod  approximately  one  mm  in  diameter  heated  electri- 
cally through  an  embedded  tungsten  wire.  He  observed  very 
slow  reaction  below  723°K,  the  melting  point  of  B2O3.  The 
reaction  rate  increased  gradually  until  vigorous  reaction 
started  near  1600-1800°K  where  the  B2O3  vaporization  rate 
became  significant.  He  Investigated  the  region  where  the 
evaporation  of  the  surface  oxides  were  rate  limiting  and 
found  the  combustion  rates  to  be  inversely  proportional  to 
pressure  and  to  obey  a temperature  dependent  Arrhenius  law 
with  an  activation  energy  of  77  Kcal/mole.  The  activation 
energy  corresponded  to  the  heat  of  vaporization  of  B2O3.  By 
adding  water  vapor  to  the  oxidant,  he  observed  an  increase 
in  the  reaction  rate  and  calculated  a change  in  the  activa- 
tion energy  to  56  Kcal/mole.  The  pressure  law  was  the  same 
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indicating  th8  reaction  was  still  controlled  by  an  evaporat- 

* 

ing  species.  This  indicated  that  if  the  surface  oxide  could 
be  removed,  the  ignition  and  combustion  characteristics 
could  be  greatly  improved.  This  would  be  true  in  the  region 
where  the  surface  oxide  vaporization  was  the  rate  controll- 
ing factor. 

Above  the  boiling  point  of  B2O3,  the  oxide  deposit 
disappeared  leaving  a clean  reacting  surface.  In  this 
region,  the  reaction  rate  was  independent  of  pressure.  Here 
boron  would  exhibit  the  ideal  heterogeneous  surface  reaction 
with  the  reaction  rate  limited  either  by  the  surface  chemical 
reaction  or  the  diffusion  of  oxygen  to  the  reacting  surface. 

Most  studies  of  particulate  boron  used  particles  or 
particle  agglomerates  from  30  to  50  /t  average  diameter. 

Cassel  (Refs  6 and  7),  Gordon  (Hef  8)  and  Macek  (Ref  9)  used 
a hydrocarbon  flame  as  the  Ignition  source.  Kelley  (Ref  10) 
and  McLain  (Ref  11)  ignited  the  particle  by  laser  energy. 
Cassel  and  Gordon  reported  on  the  gross  characteristics 
while  Macek  recorded  the  burniiig  times.  Macek  used  photo- 
graphic film  to  record  the  residence  time  of  the  particle 
in  the  flame  before  "first  light"  and  computing  the  heat 
rise  of  the  particle,  he  calculated  an  ignition  temperature 
of  near  1900°K.  Kelley  and  McLain  studied  nontreated  boron 
and  LiF  coated  particles  to  determine  whether  the  surface 
treatment  would  influence  the  chemical  klnetios.  They  re- 
ported no  observable  difference.  McLain  did  report  a factor 
of  two  Improvement  on  burning  time  for  particles  irradiated 
by  cobalt-60  gamma  radiation.  No  explanation  is  given  for 
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this  phenomenon.  The  latter  three  Investigators,  Ma&ek, 
Kelley  and  McLain,  all  report  diffusion  limited  reaction 

o 

processes  or  the  so  called  Hd-burning  law".  They  also 
assumed  an  unlimited  oxygen  environment  and  a constant 
combustion  temperature. 

Uda  and  Jones  (Befs  13  and  14)  Investigated  dense  cloud 
Ignition  In  a shock  tube.  They  found  that  ignition  and 
burning  times  of  multiple  particle  boron  systems  agreed  . 
with  the  trends  of  the  single  particle  investigations.  They 
reported  a par.ticle  size  dependence  on  Ignition  temperature 
with  0.015/c  particles  (this  particle  size  was  later  report- 
ed to  be  0.06yut  , Bef  16)  Igniting  near  1200°K  and  30-40 /t 
agglomerates  igniting  near  1900°K  at  one  atmosphere.  Gordon 
noted  the  particle  Interactive  effects  of  multiple  particles 
on  ignition  temperature.  He  reported  that  although  single 
particles  Ignited  with  difficulty  and  burned  slowly  In  a 
Bunsen  burner  type  flame,  denser  concentrations  ignited 
readily  and  burned  rapidly.  Bryant  (Bef  12)  reported  an 
Increase  in  the  time  required  to  bum  35/*-  particles  in  a 
dense  cloud  over  that  reported  by  MaSek.  However,  by  taking 
the  influence  of  oxidizer  depletion  into  account,  the 
effective  single  particle  burning  time  was  in  agreement  with 
Macek*s.  The  data  by  Maoek,  Bryant  and  McLain  are  in  agree- 
ment on  the  d^law  for  boron  particle  burning  rates  with 
approximately  40  ms  for  a 60/*.  particle  and  15  ms  for  a 
35y*  one.  This  indicates  diffusion  controlled  reaction 
rates  with  no  effects  due  to  surface  coatings  or  combustion 
environments • 
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| All  of  the  aforementioned  studies  used  ar.i  external 

i 

| means  to  Ignite  and  sustain  the  combustion  of  boron.  Bryant 

| attempted  to  establish  a steady  state,  self-sustained  flame 

s 

f but  was  unsuccessful.  He  attributed  his  lack  of  success  to 

f 

I the  large  size  particles  and  to  particle  flow  fluctuations. 

A steady  state,  self -sustained  boron/air  flame  was  success- 
fully established  by  Ortwerth  and  Wilkinson  (Ref  17).  They 

, supplied  a flow  of  Trona  amorphous  boron  particles  in  air  to 

a 1.9-centimeter  steel  combustion  tube  and  ignited  the 
boron/air  mixture  with  an  acetylene  torch.  Once  a layer  of 
glowing  oxides  and  boron  was  deposited  on  the  tube  walls  and 
the  boron  flame  was  well  established,  the  pilot  flame  could 
be  removed.  They  attributed  the  ability  to  sustain  a boron 
flame  without  the  aid  of  a pilot  flame  to  the  radiative 
energy  transfer  from  the  wall  deposits  ard  the  luminous 
flame  to  the  particles  entering  the  combustion  chamber. 

Their  system  did  not  have  a uniform  powder  flow  nor  was 
there  any  attempt  ;o  measure  the  particle  flow,  the  boron/ 
air  concentration  or  other  combustion  characteristics. 

The  theory  of  steady  state,  self -sustained  combustion  of 
solid  particles  is  not  new.  The  analysis  for  a plane  flame 
propagating  through  a dust  cloud  sustained  by  the  radiant 
heat  transfer  from  the  flame  zone  was  postulated  by  Nusselt 
(Ref  18).  This  theory  was  extended  to  include  the  effects  of 
the  heat  Conduction  from  the  fuel  particles  to  the  air  and 
the  concentration  density  of  fuel  particles  in  the  carrier 
air  by  Essenhigh  (Ref  19).  These  theories  considered  only 

an  effective  flame  temperature  heating  the  particles  to  an 

✓ 
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ignition  temperature  and  contained  no  information  on  the 
fuel  particle  temperature  history  sifter  ignition*  Because 
boron  particles  react  vigorously  at  their  self-sustained 
combustion  temperatures  and  provide  a highly  luminous  flame, 
the  steady  state  combustion  of  a dense  particle  cloud 
ignited  by  radiant  energy  from  the  reacting  particles  should 
be  possible.  However,  at  present,  there  is  no  theoretical 
analysis  for  the  temperature  history  from  the  beginning  of 
the  heat  up  zone  through  completion  of  the  reaction. 

Purpose  and  Scope 

The  purpose  of  the  present  investigation  is  to  estab- 
lish a one-dimensional  model  for  the  ignition  and  combus- 
tion of  a dense  concentration  of  boron  particles  in  air  and 
to  validate  it  experimentally.  The  theoretical  analysis  of 
the  model  is  derived  by  dividing  it  into  two  zones.  The 
first  zone  considers  the  energy  sources  and  heat  transfer 
mechanisms  to  heat  the  oncoming  gas -particle  mixture  to 
ignition  temperature.  The  second  zone  is  the  combustion 
region;  here  the  diffusion  controlled  or  d-law  chemical 
kinetic  model  for  the  single  fuel  particle,  reported  by 
previous  investigators,  was  modified  to  include  the  effects 
of  temperature  and  the  finite  oxygen  environment. 

A steady  state,  self-sustained  boron/air  flame  was 
investigated  to  evaluate  the  characteristic  parameters 
derived  in  the  analysis.  Combustion  of  both  fuel  rich  and 
fuel  lean  mixtures  were  investigated  to  evaluate  the  depend- 
ence of  particle  loading  on  the  combustion  characteristics. 
Different  combustor  tube  lengths  and  air  mass  flow  rates 
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were  used  to  evaluate  combustor  residence  time  and  flow 
characteristics  on  the  combustion*  The  different  combustor 
lengths  were  assumed  to  simulate  the  conditions  at  different 
cross  sections  of  a complete  combustor.  Trona  amorphous 
boron  was  used  as  the  fuel.  Trona  boron  is  reported  (and 
was  verified)  to  have  an  average  particle  size  of  one  micron; 
thus,  this  average  diameter  was  used  for  the  analysis. 

Data  on  oxygen  content  of  exhaust  gases  and  flame 
emissivity  taken  at  the  combustor  exit  plane  were  used  to 
deduce  the.  combustion  efficiency,  emissivity  and  temperature 
of  the  flame.  Finally,  these  data  were  compared  to  the 
combustion  and  chemical  kinetic  models  to  provide  quantita- 
tive values  for  the  characteristic  parameters  of  the  model. 

In  the  following  sections,  a one-dimensional  analytical 
model  will  be  derived  and  the  parameters  to  be  experimen- 
tally evaluated  will  be  identified.  A description  of  the 
experimental  combustion  system,  apparatus  and  instrumenta- 
tion will  then  be  presented.  The  discussion  of  the  results 
will  present  the  degree  of  success  of  this  combustion 
system  in  obtaining  high  combustion  efficiency  and 
theoretical  flame  temperatures.  These  results  will  then  be 
compared  to  the  predictions  obtained  from  the  theoretical 
model  and  the  findings  of  previous  investigators. 
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II.  THEORETICAL  CONSIDERATIONS 

There  are  many  analyses  for  radiant  heat  transfer  to 
solid  particles  suspended  in  a transparent  gas.  Simple 
radiant  heat  transfer  applies  to  a single  particle  or  a 
sparse  particle  mixture;  but  as  the  number  of  particles  is 
increased,  the  reflections  and  shielding  from  the  surround- 
ing particles  must  be  taken  into  account.  Nusselt  treated 
these  phenomena  statistically  to  obtain  an  analysis  for  a 
gas-particle  mixture  being  heated  by  radiant  heat  transfer. 

The  same  result  can  be  derived  by  assuming  an  optically  thin 

. * 
particle  cloud.  This  method  can  be  used  to  model  the  heat 

up  to  Ignition  of  a particle  cloud  but  no  analysis  has  been 

accomplished  to  characterize  the  radiant  heat  transfer  from 

a chemically  reacting  luminous  flame.  Nor  has  the  energy 

addition  due  to  radiant  wall  deposits  caused  by  wall 

reactions  been  treated  analytically. 

Combustion  Model 

The  combustion  model  for  this  investigation  was  develop- 
ed for  an  idealized  system  of  solid  fuel  particles  trans- 
ported in  suspension  by  a carrier  gas.  This  model  represent- 
ed a steady  state  gas -par  tide  mixture  flowing  through  an 
Insulated  chamber  as  depicted  in  Fig  1.  The  flow  entered 
the  chamber  at  room  temperature  at  approximately  atmospheric 
pressure.  The  suspended  fuel  particles  were  heated  to 
ignition  temperature  by  radiant  heat  transfer  from  the 
luminous  flame  and  heat  generated  by  the  reaction  of  fuel 
particles  that  adhered  to  the  chamber  walls.  During  the 
heat  up  phase,  that' is  prior  to  vigorous  combustion  of  the 
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' Figure  1.  Model  of  Combustion  System 
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particles  in  suspension,  the  carrier  gas  was  assumed  to  be 
transparent  to  the  thermal  radiation.  The  carrier  gas  was 
heated  by  convection  from  the  suspended  particles  and  the 
wall  deposits.  For  very  small  particles,  on  the  order  of 
one  micron,  the  temperature  difference  between  the  particles 
and  the  gas  has  been  found  to  be  very  small  (Ref  19 ) } so , 
for  this  model,  it  was  neglected. 

Once  the  suspended  particles  were  Ignited  to  vigorous 
combustion,  the  heat  of  reaction  on  the  surface  of  the 

particles-  provided  the  energy  source  to  heat  the  gas- 

*■  ♦ 

particle  mixture.  The  carrier  gas  could  no  longer  be  consid- 
ered transparent  to  radiation  and  radiant  energy  from  the 
fuel  particles  and  luminous  products  was  transferred  from 
the  flame  to  the  cold  particles  entering  the  chamber.  This 
provided  the  radiant  energy  source  to  heat  the  oncoming  gas- 
particle  mixture.  Thus,  the  combustion  system  can  be 
described  by  two  major  zones,  as  depicted  in  Fig  2. 

In  Zone  I,  the  heat  up  zone,  energy  was  transferred  to 
the  gas-particle  mixture  from  the  flame  zone  and  the  combus- 
tor walls;  the  fuel  particles  in  suspension  were  essentially 
nonreacting  and  the  carrier  gas  was  transparent  to  thermal 
radiation.  Zone  II,  the  flame  zone,  was  assumed  to  be 
reached  when  vigorous  combustion  was  initiated  and  the  gas- 
particle  mixture  had  become  luminous.  In  the  flame  zone, 
the  particles  in  suspension  had  ignited  and  the  energy 
release  from  this  reaction  was  the  dominant  source  of  heat 
for  the  gas -particle  mixture.  Here  the  radiant  energy  from 
the  luminous  fuel  particles  and  solid  and/or  gaseous 
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products  was  transferred  to  the  oncoming  cold  particle 
cloud. 

Assumptions 

The  assumptions  for  this  model  were  largely  standard 
simplifications  for  particle  flows  and  contain  those  stated 
by  Essenhigh  and  Csaba  (Ref  19 )•  They  were  as  follows: 

(a)  The  combustor  geometry  considered  was  a constant 
area  tube  and  hence  had  no  recirculating  zones  to  support 
the  flame;  therefore,  neglecting  turbulent  or  wall  mixing, 

the  gas -particle  mixture  was  one -dimensional  and  the  flow 

_ * 
properties 'were  constant  across  any  given  cross  section. 

(b)  The  boron  particles  were  spherical  of  uniform  size 
with  the  number  of  particles  n in  a given  gas  volume  a 
constant. 

(c)  The  gas  and  particle  temperatures  and  velocities 
were  identical. 

(d)  The  volume  occupied  by  the  solid  particles  was 
negligible  compared  to  the  gas  volume;  thus,  the  system 
velocity  was  considered  to  be  that  of  the  gas  without 
particles. 

(e) i  The  gas  properties  did  not  vary  significantly 
from  those  of  air. 

(f)  The  viscous  drag  on  the  combustor  wall  was 
neglected. 

(g)  The  thermal  conductivities  and  specific  heats  of 
the  gases  and  particles  were  constant  over  the  temperature 
rar.ges  concerned. 

\ 
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(h)  In  Zone  I,  the  solid  particles  within  the  mixture 
received  energy  from  the  flame.  The  flame  was  assumed  to  be 
at  an  effective  temperature  Tp. 

(i)  The  combustor  tubes  were  constructed  of  high 
thermal  resistant  materials;  thus,  the  inner  wall  of  the 
combustor  was  considered  to  be  surrounded  by  an  adiabatic 
Jacket  and  all  the  energy  generated  by  the  heat  of  reaction 
on  the  combustor  walls  during  heat  up  to  ignition 
temperature  was  returned  to  the  flow  system. 

(J)  .The  mass  of  fuel  and  air  that  burned  on  the  chamber 

♦ « 
walls  was  Jre  turned  to  the  flow  as  products. 

(k)  Before  reaction  of  the  suspended  boron  particles 
became  an  important  source  of  energy,  the  air  surrounding 
the  particles  was  transparent  to  thermal  radiation.  It 
received  its  energy  by  convection  from  the  particles. 

(l)  The  gas-particle  mixture  behaved  like  an  optically 
thin  absorbing-emitting  gray  gas  in  Zone  I. 

(m)  In  the  boron/air  reaction,  only  boron  and  oxygen 
took  part  in  the  reaction  with  B20^  as  the  primary  product. 
Definitions 

Elemental  boron  is  difficult  to  obtain  in  a completely 
pure  state  and  the  boron  used  in  this  experimental  investiga- 
tion, Trona  amorphous  boron,  has  approximately  10#  inert 
impurities  by  weight.  Table  III  of  Appendix  J.  There  are 
times  when  the  value  for  the  mass  of  boron  will  be  used  in 
the  analysis  instead  of  that  of  the  fuel.  The  subscript  f 
for  fuel  will  be  used  when  the  impurities  in  the  fuel  must 
be  considered  to  account  for  the  total  mass  of  the  fuel; 
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when  only  the  mass  of  boron  is  to  be  considered,  the  sub- 
script B will  be  used.  For  a fuel  that  has  a fuel-to-boron 
mass  ratio  mf/mB  = 1 +A»  the  fuel  is  said  to  have 
A/(l  +A  ) x 100#  impurities.  Examples  of  where  this  will 
arise  are  in  the  definition  of  combustion  efficiency  >)  and 
fuel-to-alr  mass  ratios.  From  Appendix  A,  the  definition  of 
the  combustion  efficiency  V}  is 


yi  = M.- 


(A-6) 


But  A rag  =*  Amf  since  boron  is  the  only  active  element  in 
the  fuel,  therefore 


T)  - Any  _ (W f).  Atrif  _ f/+A)AW f 

(m„\  (me\  (Wf).  (mf>. 


(A-6a) 


Fuel-to-air  ratios  are  defined  as 


f/fK-  (JZk).  = (t+A)(Jh£*  = 0 + A)(%) 


For  the  stoichiometric  condition 


07a)*  = (i+a)(B/a)s 


The  F/A  and  B/A  ratios  are  interchangeable  in  the  definition 
of  the  equivalence  ratio  <f>  . 


^ = ClMd  CM!  = (M) 

IW>S  (l+A)  (0/A)s  ( F/A), 
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The  completeness  of  the  boron/oxygen  reaction  can  be 
measured  by  the  depletion  of  the  oxygen  as  well  as  the  boron. 
For  fuel  lean  combustion  { ^ < 1),  the  boron  would  be 
completely  consumed  when  the  reaction  was  completed  and  the 
combustion  efficiency  would  range  from  zero  to  one.  For 
fuel  rich  combustion,  the  oxygen  would  be  depleted  when  the 
reaction  was  complete  but  unconsumed  fuel  would  remain. 

Thus,  the  combustion  efficiency  would  always  be  less  than 
unity,  ^jaax  = 1/^  . . However,  in  this  case,  the  oxygen 
depletion. parameter  'p  defined  in  Appendix  A as 

f ynoj/fao,)'  *£2s  (A-8) 

would  range  from  zero  to  one.  Equation  A-8  defined  the 
oxygen  depletion  parameter  and  indicated  the  relation  be- 
tween the  combustion  efficiency  and  the  oxygen  depletion 
parameter  to  be 


syj/  s £ Yj  (A— 8) 

It  -is  seen  from  Eqs  1,  2 and  3 that  the  equivalence 
ratio  <f>  is  defined  from  the  Initial  conditions  of  the 
gas-particle  mixture  and  is  therefore  a constant.  However, 
the  combustion  efficiency  and  oxygen  depletion  parameters 
define  the  completeness  of  the  combustion  process  and  are 
functions  of  time  or  combustor  station  for  a steady  process. 


DS/ME/74-2 
Control  Volume 


The  elemental  control  volume  of  the  combustor  model. 

Pig  3*  was  used  to  apply  the  conservation  equations  of  the 

flow  and  to  determine  the  range  of  the  variables  encountered 

In  this  combustion  system.  The  Initial  conditions  of  the 

flow  Included  the  temperature  to  be  at  room  temperature  and 

the  velocity  In  the  low  subsonic  range,  while  the  temperature 

condition  at  the  combustor  exit  was  at  or  greater  than  the 

ignition  temperature.  The  ignition  temperature  was  reported 

by  all  investigators  to  be  near  1900°K.  For  complete 
* ♦ 
combustion' however,  the  exit  temperature  was  expected  to  be 

near  the  adiabatic  flame  temperatures  shown  on  Fig  18. 

Conservation  Equations 

The  conservation  equations  for  a gas-particle  mixture 
were  developed  In  Appendix  C.  The  assumptions  for  these 
equations  included  one-dlmenslonal  flow,  no  viscous  forces, 
constant  area  duct  and  Identical  velocity  and  temperature 
for  all  phase  states.  Since  the  velocity  of  all  phase 
states  was  the  same,  the  mass  within  an  elementary  control 
volume  at  distance  x from  the  combustor  entrance  was  propor- 
tional to  the  mass  in  an  identical  volume  at  the  combustor 
entrance.  Therefore,  even  when  combustion  takes  place  and 
the  chemical  species  or  phase  states  change,  the  mass  of  the 
elementary  (atomic)  species  remains  constant  within  any 
given  volume. 
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The  term  LE/g.  Cp  describes  the  volumetric  specific  I 

Hi  I"  ”f  *’ft 

heat  of  the  multiphase  flow.  This  term  was  evaluated  in 
Appendix  D and  was  determined  to  be 

where  © s / + 0. 3/  <j>  (D-8) 

The  pressure  term  in  Eq  C-21  is  approximately  equal  to  the 
change  in  the  kinetic  energy.  For  flow  velocities  less  than 

. - ..d 
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sonic,  the  change  due  to  kinetic  energy  is  small  when  compar- 
ed to  the  combustion  heating  and  heat  transfer  from  the  flame; 
thus,  it  can  be  neglected.  The  heat  transfer  term  consists 
of  both  radiant  and  conductive  energy  transfers  with  both 
longitudinal  and  radial  components • The  radial  energy 
transfer  from  the  combustor  walls  was  considered  as  a single 
energy  source,  neglecting  angular  dependence  Eq  C-21  can  be 
written 


Cp  elXT 

a * *"•  d*  ~ a,r  “ tx 


7 Jr  ^ ^ + ^ 


Axial  Heat  Transfer 

Even  though  the  axial  transfer  of  energy  is  a combined 
process  of  conduction  and  radiation,  one  of  these  processes 
may  be  much  larger,  Just  as  a significant  energy  change  by 
heat  transfer  overshadows  the  changes  due  to  kinetic  energy. 
To  determine  the  relative  importance  of  the  axial  conduction 
versus  the  radiation  heat  transfer,  the  surface  integral 
describing  heat  transfer  in  Eq  C-14  can  be  evaluated  over 
the  volume  encompassing  the  entire  heat  up  zone.  The 
boundary  conditions  for  this  volume  are  ambient  temperature 
TQ  at  the  combustor  entrance  and  a temperature  at  or  greater 
than  the  self -sustained  ignition  temperature  T^  at  the 
boundary  between  the  heat  up  zone  and  the  combustion  zone. 
Past  Investigators  have  established  T^  to  be  near  1900°K. 

The  radiant  heat  transfer  at  these  surfaces  would  be  <f  T0^ 
and  Tp^„  respectively,  where  the  combustion  temperature  Tp 
is  an  average  temperature  of  the  flame  zone.  Integrating 
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the  heat  transfer  surface  Integral  between  these  extremes 

» 

and  noting  that  the  surface  areas  are  constant  yields 


-JJr^-J/i-js-f/%-js 

=-  [etc  ^+tysl]  - [^4+<jf4] 


Taking  the  ratio  of  the  radiant  transfers  to  conductive 
transfers  In  Eq  5 and  assuming  = Tp  = 2000°K  with 
O’  = 1.355  * 10"12  cal/cm2-sec-K^and  J^lr  = 1.2  x 10",*f 
oal/om-sec-K 


W : o-  (Tx4-TS)  (Xx-K  ) 

/ p r-r  _ -r  ) 

•air  ' U '9  ' 

A/ 

0i)(/.7)  (/O-'l 


- (x*-X')xio* 


When  the  heat  up  distance  (x^  - xQ)  is  on  the  order  of  a 
centimeter,  the  radiant  transfer  is  much  greater  than  that 
of  conduction.  To  obtain  a self -sustained  flame,  Ortwerth 
and  Wilkinson  (Ref  17)  required  a minimum  combustor  length 
of  15  centimeters,  much  greater  than  one  centimeter.  There- 
fore, It  was  concluded  that  the  conductive  heat  transfer 
could  be  neglected. 

The  radiant  energy  transfer  into  en  elemental  control 
volume  of  the  gas-particle  mixture,  in  which  Eq  k applies, 
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is  absorbed  or  reflected  by  the  solid  particles.  The 
particles  also  emit  radiant  energy  which  is  reflected  by  the 
surrounding  particles.  In  the  combustion  zone,  the  solid  or 
gaseous  combustion  products  could  also  contribute  to  the 
absorptions,  emissions  and  scattering  of  the  radiant  energy. 
Thus,  whether  the  elementary  volume  is  within  the  heat  up  or 
combustion  zone,  the  gas-particle  mixture  will  behave  like  an 
absorbing,  emitting  and  scattering  media.  Applying  the  above 
boundary  conditions,  the  axial  radiant  energy  transfer  and 

its  derivative • were  developed  in  Appendix  E for  the  gray  gas 

«•  « 

media.  Invoking  the  optically  thin  approximation  assumption, 
-dqr/dx  can  be  written 


^ zA*  [r?+l*-2T*] 


(E-9) 


Where  Tp  is  the  effective  temperature  of  the  flame  zone,  TQ 
the  Initial  temperature  of  the  gas -particle  mixture  and  O' 
the  Stefan-Boltzmann  constant.  While  k and  Ji  are  defined 
as  the  volumetric  absorption  coefficient  end  the  mass 
absorption  coefficient,  respectively. 

Absorption  Coefficients 

The  combustion  zone,  Zone  II,  was  defined  as  the  region 
where  the  cloud  reactions  became  important  and  the  brilliant 
green-white  flame,  that  has  been  described  by  all  investi- 
gators of  boron  combustion  appeared.  Kasken  and  Millikan 
(Ref  23),  John  (Ref  24)  and  Tisher  and  Scheller  (Ref  25) 
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reported  the  luminosity  of  the  flame  to  eminate  from  the 
gaseous  reaction  products.  John  identified  the  principle 
visible  emitter  as  B02^  Thus,  while  the  gases  in  Zone  I 
were  transparent  to  radiation  and  the  particles  made  up  the 
absorption  coefficient,  the  gaseous  products  in  Zone  II 
added  to  the  absorption  coefficient.  This  lead  to  the  model 
for  the  volumetric  absorption  coefficients. 

For  Zone  I, 

..  4*4. 


Where  k^  was  .the  absorption  coefficient  due  to  the  particles 
* ■% 
in  the  gas -particle  mixture. 

For  Zone  II, 

A - Ap  + -Afr 

Where  kg  ' ns  the  reaction  products  contribution  to  the  total 
mixture  absorptivity  and  must  be  modeled  with  the  assistance 
of  experimental  results  from  this  investigation. 

The  Particle  Volumetric  Absorption  Coefficient 

The  volumetric  absorption  coefficient  of  the  suspended 
particles  was  modeled  after  the  reciprocal  of  the  mean  free 
path  of  a gas  from  molecular  theory  and  was  proportional  to 
the  total  cross-sectional  area  presented  by  the  particles 
per  unit  volume. 


: 


■A 


$ 


Ap  = Z Y]  TTc/ 
4 


(6) 


Where  d is  the  average  diameter  of  the  fuel  particles,  n the 
number  density  of  particles  per  unit  volume  and  Z is  a 
dimensionless  proportionality  factor.  The  number  n can  be 
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calculated  by  dividing  the  mass  of  the  solid  fuel  per  unit 
volume  m^.  by  the  mass  per  particle.  For  spherical  particles 
this  is 


SgTT^3 


(7) 


Where  &■%  is  the  density  of  solid  boron  and  dQ  the  initial 
particle  diameter.  Equation  7 can  be  described  in  terms  of 
the  density  of  the  carrier  air  ^air  and  the  equivalence 
ratio  4*  as 


/£„  j (JMs  t 

Sq,  TT  4* 


(8) 


Substituting  Eq  8 into  Eq  6,  the  volumetric  absorption 
coefficient  of  the  particles  is  then 


ip  = *2  /g-  V/A\  4.  (9) 

2 $3  Jo 

The  combustion  efficiency  Y\  was  de**ned  in  Eq  A-6a  as 

(1  +A)(mf.  -mf.)/m.P  . For  spherical  particles,  it  can  be 
■*•0  1 Ao 

described  in  terms  of  the  particle  diameter  ratio 


Thus,  the  volumetric  absorption  coefficient,  Eq  9,  may  be 
written  as 

4 = % p„  2 <m  * [ / - % . 

&B  Jo 
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and  the  mass  absorption  coefficient  as 

*p=  % i (r/A)*  * (l~n')% 

$8  d0 

Introducing  Y)'=  >)/{l  +A). 

Using  the  values  dQ  = 1 x lO^cm,  &%  - 2.34  gm/cc, 
(P/A)s  = 0.115  and  £f=  1.355  x 10“12  cal/cn^-sec-K**',  the 
coefficient  of  the  radiant  heat  transfer  due  to  the  particles 
may  be  evaluated  as  . 

. . iXfO-.  zzxis'1 4>  u-y)Vi  (a) 


Or  expressed  in  terms  of  the  oxygen  depletion  parameter  Y' 
from  Eq  A-8,  Eq  11  may  be  written 


% 


22  X/09<t>(/-%')  cJ/f„-k4-sec  (12) 


Where  (f>  = (1  +A)<f>. 

Chemical  Heat  of  Reaction 

The  energy  addition  per  unit  volume  me  to  the  chemical 
energy  liberated  by  the  combustion  process  Qchem,  in  Eq  4, 
was  equal  to  the  heat  of  formation  of  the  products  per  unit 
mass  of  fuel  multiplied  by  the  rate  of  fuel  consumption,  thus, 





4*m  = -<  P 


= (”klA4$ 
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Where  (mB)0  was  the  mass  of  boron  atomic  species  in  a unit 
volume  and  ^Jh£  the  energy  liberated  by  the  production  of 
B2O3  per  unit  mass  of  boron. 

Equation  13  can  be  described  in  terms  of  the  density 
of  the  carrier  air  since  the  density  of  the  air  is  the  mass 
of  air  per  unit  volume.  Thus,  the  mass  of  boron  in  a unit 
volume  divided  by  the  air  density  is  the  boron/air 

ratio  (B/A) 

.v.. 

• ' • . . • ■»  . • ■ , ..  • • 
. • ...  , r . ' • 

Expressed  in  terms  of  the  stoichiometric  ratio  and  the 
equivalence  ratio  Eq  13  may  be  written  as 

Q'hen Zr  (e/A%4>^(  jl  (14) 

The  measurement  of  the  oxygen  content  of  the  exhaust 
gases  can  be  accomplished  with  ease.  However,  the  deter- 
mination of  the  percentage  of  boron  that  has  reacted  would 
be  much  more  difficult.  Therefore,  Eq  14-  can  be  written  in 
terms  of  the  oxygen  depletion  parameter.  The  equivalence 
ratio  can  be  combined  with  the  differential  of  the  combus- 
tion efficiency  in  Eq  14  and  the  chemical  energy  release 
rate  may  be  written  as 

Qchem  = sf,r  (B/a)s*A 
s = /g,  (m)s^f 

/ 
s 
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In  Zone  I,  the  heat  up  zone,  the  radial  energy  addition 
from  the  combustor  walls  must  be  considered.  The  source  of 
this  energy  was  chemical  reaction  from  particles  reacting  on 
the  combustor  walls.  But  since  the  reactions  took  place  on 
the  walls  outside  of  the  control  volume,  the  energy  must  be 
transferred  back  into  the  control  volume.  This  would  be 
accomplished  by  a complicated  process  of  the  combination  of 
radiant,  convective  and  mass  transfers.  To  simplify  the 
mathematical  relations  for  the  process,  the  radial  energy 
flux  generated,  by  the  wall  reactions  was  considered  as  part 
of  the  volumetric  chemical  heating.  The  energy  generated  by 
the  chemical  reaction  in  Eq  4 can  then  be  written  as  the  rate 
of  chemical  energy  released  by  the  combustor  wall  reaction 
plus  the  rate  of  the  chemical  energy  released  from  the 
particles  suspended  in  the  air  within  a control  volume. 

r Qh em^w*  Qch ey»')c/oltJ 

This  relation  may  be  stated  in  terms  of  the  oxygen 
depletion  rate;  therefore,  Eq  15 » which  is  a general 
expression,  can  be  written  as 

Q<U,  - (B/a\  + 

Where  (d^/dt)  was  the  rate  of  oxygen  depletion  on  the 

W 

combustor  walls  and  (d^/dt)olou(i  the  oxygen  depletion  rate 
for  the  suspended  particles  reacting  in  the  carrier  gas. 

\ 
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The  heat  of  formation  for  the  product  BgO-j  decreases  as 
the  temperature  Increases  due  to  the  change  In  the  phase  of 
BgO-j  from  liquid  to  gaseous.  To  reflect  this  change  In  the 
heat  of  formation  term,  Ah°  was  expressed  as  a function  of 
the  heats  of  formation  for  the  liquid  and  gaseous  products. 


A-flf*  A A ^£rzo3lA^-f^)  A 


(17) 


was  the  partial  pressure  of  B20o  in  atmospheres 

b2°3 


Where  P. 

at  ambient,  pressure.  The  value  for  Pv 
Edelman  and  Economos  (Ref  21)  was 


B2O3 


postulated  by 


r„  = e 


*P 


[ \q-  44coo/t] 


The  values  for  Ah°(j{)  and  Ah^(g)  were  taken  from  the 
JANAF  tables  (Ref  22): 

4<tf)  = 13.1  Kc«*/j„b 

AJlfCtf)  - K’po^jrig 

thus 

A*(=  /3.7  - Exp  [fl-Hoofr]  4.T4  HutA^  (18) 
Wall  Reaction  Rate  Model 

The  reaction  on  the  combustor  walls,  (dl^/dt)^,  would 
be  proportional  to  the  rate  at  which  the  particles  impact  the 
walls,  the  tempera turo  of  the  combustor  walls,  the  temperature 
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of  the  particles  before  impact  and  the  condition  and  thick- 
ness of  the  oxide  deposits  on  the  walls.  It  could  also  be 
affected  by  the  velocity  of  the  air  across  the  deposit 
surface.  The  temperature  of  the  walls  and  the  condition  of 
the  wall  deposits  would  be  difficult  to  establish  and  the 
velocity  of  the  air  will  increase  as  the  particle  tempera- 
ture increases.  Therefore,  for  a simplified  relation,  the 
wall  reaction  rate  was  assumed  to  be  proportional  to  the 

equivalence  ratio  £ and  the  particle  temperature, 

• • * 

(#1  = ?T*  (19) 

Where  g is  a constant  of  proportionality  with  units  of 
(sec-K)"1.  In  Zone  II,  the  energy  generated  by  the 
reactions  of  the  particles  suspended  in  the  gas-particle 
mixture  was  assumed  to  be  the  dominant  chemical  energy 
source  and  the  wall  reactions  were  neglected. 

Particle  Reaction  Rate  Model 

The  reaction  rate  for  the  wall  combustion  processes 
was  modeled  from  experimental  results  while  the  cloud 
combustion  process  for  the  suspended  particles  was  modeled 
from  single  particle  combustion  models  of  previous  investi- 
gators (Hefs  9,  10,  11  and  12).  The  present  model  for 
suspended  particles  differs  from  past  models  in  that  the 
reaction  temperature  was  assumed  to  Increase  as  the  combus- 
tion process  proceeded.  All  other  models  have  assumed  a 
constant  reaction  temperature.  The  constant  temperature 
assumption  may  be  correct  for  the  nonadiabatlc  combustion  of 
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particles  radiating  their  thermal  energy  to  the  surroundings 
at  ambient  temperature  or  for  liquid  particles  burning  with 
a diffusion  flame  where  most  of  the  thermal  energy  is 
required  to  vaporize  the  liquid  fuel.  However,  the  present 
combustion  system  was  a dense  cloud  of  solid  particles  burn- 
ing in  an  insulated  chamber  to  approximate  an  adiabatic 
combustion  process.  The  boron  particles  may  be  molten  at 
the  combustion  temperatures  but  the  volatility  of  the 
particles  was  assumed  to  be  so  low  that  the  combustion 
process  may  be-  represented  as  a surface  reaction  as  depicted 
in  Fig  4. 

For  a solid  fuel  reacting  with  a gaseous  oxidizer,  the 
rate  of  consumption  of  the  solid  material  is  proportional 
to  the  surface  area  of  the  solid  and  the  concentration  of 
the  oxidizer  in  the  surroundings.  For  a spherical  solid 
particle  burning  in  air,  the  empirical  first-order  diffusion 
limited,  reaction  rate  is  described  in  Appendix  I as 


JhflB  _ 4 yj  7T  Y3-  Co  K(T 2 

dir  “ t 


(20) 


Where  m^  was  the  mass  of  boron  per  unit  volume,  n was  the 
number  of  fuel  particles  per  unit  volume,  r the  radius  of 
the  particle,  C0  the  concentration  of  oxygen  in  air,  Mg  the 
molecular  weight  of  boron,  and  K(T)  a temperature  dependent 
proportionality  coefficient  in  moles/cm-sec. 

The  mass  of  a spherical  particle  of  boron  can  be 
expressed  in  terms  of  the  particle  diameter  d,  thus,  Eq  20 
may  be  written  as  ' 
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Figure  4 
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Diffusion  of  Product 
Oxides  away  from  the 
Particle  Surface 


- — Diffusion  of  Oxygen  to 
the  Particle  Surface 


Heterogeneous  Surface  Reaction  Model 
for  Solid  Shperical  Particles 
Suspended  in  the  Carrier  Gas 
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= ZftTTdc^rte  Krt 


(21) 


For  constant  temperature  reaction  and  an  unlimited  source  of 
oxygen,  Eq  21  can  be  integrated  to  obtain 

<£-/  = t_Coi  Mn  Kffl  ft-t-l 

<^B 


(22) 


Where  y5dlff  = B(CQ  )Q  Mg  K(T)/£fi.  Equation  22  is  the  d-law 

for  diffusion  limited  reaction  of  a fuel  particle  or  droplet 

and  /^iff  Is  defined  as  the  d-law  evaporation  rate  constant. 

For  gas-particle  mixtures  near  the  stoichiometric 

mixture,  the  oxygen  concentration  Is  no  longer  constant  but 

varies  with  the  combustion  efficiency.  Nondimens lonalizing 

Eq  21  by  dividing  this  relation  by  dQ^  and  multiplying  and 

dividing  by  the  Initial  concentration  of  oxygen  (CQ  )Q,  one 

2 

obtains 


--  w°kfr,^)[§;J  («) 

For  a constant  combustion  temperature,  only  the  oxygen 
concentration  and  the  particle  diameter  ratios  change  with 
time.  These  ratios  can  be  expressed  as  functions  of  the 
combustion  efficiency  by  using  the  relations  developed 
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i:.  Eqs  10  and  A-8.  Substituting  these  relations  into  Eq  23i 
it  can  be  expressed  In  terms  of  the  combustion  efficiency  as 

jl ) . |2  a )rt,  (<«.)„  fr(T>  (l-rfY*  (t- 4*  ) (24) 

<tt  ' Sb  ct1 

Now,  to  extend  this  analysis  to  a temperature  dependent 
model,  the  proportionality  coefficient  K(T)  was  expressed 
as  an  Arrhenius  equation 


KM  = A e“6l'/RT  (25) 

Where  Ea  and  A are  the  Arrhenius  activation  energy  and  the 
Arrhenius  pre-exponential  constant,  respectively.  The  pre- 
exponentia1  constant  A has  a temperature  dependence  in  the 
more  exact  theories;  but  unless  very  accurate  experimen- 
tal measurements  are  made,  this  is  difficult  to  detect.  The 
temperature  dependence  of  the  exponential  term  normally 
dominates  in  most  reaction  systems;  therefore,  the  term  A 
was  considered  as  a constant. 

Replacing  the  proportionality  factor  K(T)  in  Eq  25  with 
the  Arrhenius  expression,  the  temperature  dependent  particle 
reaction  model  for  a diffusion  limited  reaction  rate  and  a 
finite  oxygen  source  expressed  in  terms  of  the  combustion 
temperature  and  combustion  efficiency  was 

_£>  - I2(  I* A) (C„r\  H J (/-*!')  M E*p[-  E*/f?T]  (26) 

\ 
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Equation  2 6 may  ’ <tten  in  terms  of  V as 

JV  _ iZtfoJ.HaA  <V  (l-V/pfiCl-V)  E*p[-E./RTj  (27) 
J<  ’ 4* 

Where  <f/ (1  +A)  4 and  (C0  )0  was  the  initial  oxygen 

2 

concentration,  moles  02/moles  air,  equal  to  0.21.  This 
relation  was  for  the  particles  reacting  in  the  particle 
cloud  with  the  reaction  limited  by  a diffusion  process. 
Equation  27  can  be  simplified  as 

drf  _ ft'  tfCl-V/p)  3 (l.-V) E*p  ~ “)]  (28) 

di  " 

Where  A»dlff  = 12(C0  )0  Mfi  A Exp  [-Ea/R232o]/  SBdQ2  with 
units  of  (sec)"1. 

The  Residence  Time  Parameter 

The  reaction  rate  model  for  the  combustion  process, 

Eqs  19  and  28,  has  been  defined  as  a function  of  residence 
time  t;  while  the  Increase  in  temperature  of  the  combustion 
process,  Eq  4,  was  expressed  as  a function  of  the  axial 
distance  from  the  combustor  entrance.  For  a steady  process, 
the  relation  between  the  axial  distance  and  time  is 

dX=  U di 

or 

JZ/Uo  “ (vAO  <Li  ( 29 ) 

Using  the  continuity  equation,  Eq  29  can  be  written 
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(^/4tW£  - dj£  - J%  =:  ol  % (30) 

^ c&n 

Where  0 ~ £(/&)*  *‘/0%„'VC)o  ] is  defined  as  the  residence  time 
parameter. 

The  velocity  ratio  can  be  written  as  a function  of 
temperature.  From  the  analysis  In  Appendix  F 


/*  ( /+fr/A')y  fa  JL  (F-12) 

/ + U+F/A^mlX  J T° 

If  the  Mach  number  Is  neglected,  significant  error  would 
not  be  introduced  for  Mach  numbers  much  less  than  one. 

Thus,  for  the  present  Investigation,  Eq  30  becomes 


%.  it  S X Jt  = (3D 

V.  T„a 

The  energy  equation  Eq  4,  along  with  the  relations 
developed  for  the  radiant  heat  transfer  Eq  E-9  and  the 
chemical  energy  release  Eq  16,  comprises  the  one-dimensional 
combustion  model  for  the  gas -particle  mixture.  These 
equations  apply  from  the  beginning  of  the  heat  up  zone 
through  the  flame  zone.  This  model  may  be  written  in  terms 
of  the  residence  time  parameter  g as 


G/?,rcP  JJ.  - 2/°ir  X * ( T*"'‘  T*  “ 

e,r  rair 


% 


\ 


(32) 
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Here  we  have  neglected  the  conductive  heat  flux.  The  first 
term  on  the  right  represents  the  gradient  of  the  axial 
radiant  heat  flux,  the  first  portion  of  the  second  term 
represents  the  radial  heat  flux  and  the  last  portion  repre- 
sents the  volumetric  heating  due  to  the  combustion  of  fuel 
particles  in  suspension. 

Multiplying  Eq  32  by  TQ/T  (noting  that  u/uQ  = T/TQ)  and 

dividing  by  9A,  c_  , the  combustion  model  for  the 

*alr 

temperature  rise  as  a function  of  the  residence  time 
parameter  is  . 

il  . 2_XrZ  z*-zr>h 

d?  &Cp  T ©Co 

'air  loir 


The  Combustion  Model  for  Zone  I 

In  Zone  I,  it  was  assumed  that  the  gases  were  trans- 
parent to  the  radiant  heat  transfer  and  only  the  solid 
particles  in  the  mixture  absorbed  the  radiation.  The 
absorption  coefficient  X Is  only  that  due  to  the  particles. 
The  dominant  chemical  reaction  in  this  zone  was  the  reaction 
at  the  combustor  wall  while  the  reactions  on  the  suspended 
particles  became  increasingly  important  as  the  temperature 
approached  the  self-sustained  ignition  temperature.  The 
implicit  functions  and  constants  of  Eq  33  can  now  be 
evaluated  and  expressed  as  functions  of  the  equivalence 
ratio,  the  oxygen  depletion  parameter  and  temperature. 

Using  the  numerical  values  for  c and  (B/A)_  of  0.257 

pair  s 


, qrgec  "vtTr;^>£  rx?  ^ ^^^r<~  \^v-fz  ■*<%  v A^T-  wr;--4','.'  ^ ?vY  ^ ??r- 

'*46  "*  ~ ' ' 
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cal/gm-K  and  0.1046,  respectively,  and  the  expression 
2d(p^  developed  In  Eq  12,  then 

ZXgtrZ  - g.  33Z  //o"i’(z5  (z-V/fif*  (sec -k1)  1 (34) 

%, 

Using  the  numerical  expression  for  the  chemical  heat  of 
formation  for  BgO-j,  Eq  18,  the  coefficient  for  the  chemical 
energy  term  in  Eq  33  was  calculated  to  be 


(B/a\  . = ftso  - Exp  [ i1-4+*°*/t\/8*o  °H  (35) 

Co  * 
ratr 

The  simplified  expression  for  the  chemical  reaction 
rate  on  the  wall  was  developed  from  the  relation  established 
In  Eq  19  and  can  be  written  in  terms  of  the  residence  time 
parameter  £ as 


Jk(H) 

T t dtl 


(36) 


The  model  for  the  reaction  rate  of  the  particles 
suspended  In  the  carrier  air  Is  described  in  Eq  28.  This 
relation  can  be  written  in  terms  of  the  residence  time 
parameter  £ as 

(*£)"  j=  **&£(!- v/# f3( ' I -V)Ev [<M&»T  7)] < 3? ) 

Equation  33  is  the  general  expression  of  the  model  for 
the  combustion  of  a gas -particle  mixture.  This  relation, 
together  with  Eqs  34  through  37 » represents  the  gas -particle 
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mixture  before  vigorous  combustion  takes  place.  This  model 


accounts  for  reaction  of  fuel  particles  on  the  combustor 


walls,  the  initial  stages  of  the  energy  release  from  the 


particles  in  suspension  and  neglects  axial  conductive  heat 


transfer.  The  carrier  air  is  assumed  transparent  to 


radiation  by  considering  absorption  and  emission  of  radiant 


energy  by  the  fuel  particles  alone. 


The  Combustion  Model  for  Zone  II 


In  Zone  II,  the  suspended  particle  reaction 


(dV'/dl' ) cloud  was  fche  PrinoiPal  reaction;  the  wall  reaction 
(df/dZ  )w  was  assumed  to  play  a minor  role  in  the  reaction 


processes  and  was  neglected.  The  luminous  products  added 


to  the  luminosity  of  the  flame  and  were  included  in  the 


absorptivity  term  of  the  heat  transfer.  A control  volume 


within  this  zone  was  assumed  to  be  in  thermal  equilibrium 


with  the  rest  of  the  flame;  its  only  thermal  exchange  was 


with  the  cold  gas-particle  mixture  approaching  the  flame. 


For  Zone  II,  the  equation  for  the  rate  of  temperature  rise 


Eq  33  was  then 


JT  <r  (Af  ) . ,,Q, 

7i- 


Where  the  relations  established  in  Eqs  34,  35  and  3?  are 


still  valid. 


The  parameters  affecting  the  combustion  model  have  been 


identified  and  its  theoretical  analysis  may  be  accomplished 


with  Eqs  33  through  38.  Next,  an  experimental  parametric 


study  of  the  steady  state,  self -sustained  boron/air 
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combustion  will  be  required  to  determine  the  characteristics 
of  the  boron  combustion  process.  The  parameters  to  be 
evaluated  ares  the  fuel/air  equivalence  ratio  (£> , oxygen 
depletion  parameter  /y/,  the  local  gas-particle  temperature 
T,  the  Arrhenius  equation  activation  energy  Ea  and  pre- 
exponential constant  A,  the  effective  flame  temperature  Tp, 
the  combustor  wall  reaction  rate  constant  g.  the  particle 
absorption  proportionality  factor  Z,  the  combustion  product 
contribution  to  the  flame  absorptivity  v v and  the  combus- 

o 

tor  residence  time  parameter  % . 

The  parametric  study  was  accomplished  at  different 
values  of  the  residence  time  parameter  % and  the  equiva- 
lence ratio  <j>  • The  residence  Lime  parameter  f;  was  varied 
by  changing  the  combustor  length  L and  the  air  mass  flux. 

The  equivalence  ratio  <f>  was  determined  by  metering  the  fuel 
flow  for  the  corresponding  air  mass  flow.  The  oxygen  deple- 
tion parameter  and  the  gas-particle  temperature  were 
determined  for  the  various  values  of  % and  <f> . The  oxygen 
depletion  values  were  measured  from  exhaust  gas  samples  by 
mass  spectrography  and  the  temperatures  were  evaluated  from 
absorption-emission  data.  The  rest  of  the  parameters  were 
evaluated  using  the  combustion  model  to  determine  the  best 
correlation  to  the  measured  data. 


\ 


39 


■jt.1  . 


DS/ME/74-2 

III.  APPROACH  AND  EXPERIMENTAL  APPARATUS 
A system  was  developed  that  permitted  a high  degree 
of  control  over  the  fuel  particle  flow  for  a self-sustained 
boron/air  combustion  flame.  The  system  consisted  of  dry, 
filtered,  compressed  air  and  boron  particles.  The  air  supply 
was  metered  through  a central  duct.  The  air  was  then  divided 
with  a portion  of  the  air  directed  through  a particle  mill 
with  the  main  flow  of  air  going  directly  to  a combustion 
chamber.  The  air  directed  through  the  particle  mill  re- 
joined the.  main  flow  upstream  of  the  combustion  chamber. 

Boron  fuel'  was  prepared  by  pressing  Trona  amorphous  boron 
into  cylindrical  pellets.  The  pellets  were  placed  into  a 
fuel  cylinder  and  fed  into  the  particle  mill  at  a given  rate 
to  obtain  the  desired  flow  rate  of  boron  particles.  The 
particles  were  entrained  in  the  air  flow  that  was  directed 
through  the  particle  mill  and  carried  to  the  main  flow.  The 
particle  laden  flow  was  injected  longitudinally  along  the 
center  line  of  the  inlet  duct  46  centimeters  upstream  from 
the  combustion  chamber  to  allow  for  mixing  the  particles  and 
air  before  they  entered  the  combustor.  The  boron/air  mixture 
was  initially  ignited  by  a hydrogen/air  pilot  flame.  Once 
the  boron  combustion  was  established,  the  pilot  flame  was 
turned  off  and  the  self -sustained  boron  flame  would  stabilize 
at  near  atmospheric  pressure.  After  the  flame  had  reached 
steady  state  combustion,  data  on  flame  tempera ture  and 
exhaust  gas  oxygen  content  were  taken.  A diagram  of  the 
test  apparatus  is  shown  in  Fig  5* 


, tv  — ' 
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Test  Matrix 


The  characteristics  of  self -sustained  combustion  at 


near  atmospheric  pressure  were  investigated.  The  combustion 
chambers  used  were  constant  area  tubes  with  no  restrictions 
at  the  exit.  Combustor  tubes  up  to  76  centimeters  in  length 
were  used  to  vary  the  combustion  residence  time.  The  length 
of  the  combustor  was  measured  from  the  exit  plane  back  to 
the  igniter  tube  tip.  The  fuol-to-air  loading  effects  were 
determined  by  varying  the  fuel/air  equivalence  ratio  from 
0.5  to  2.0.  The  effects  due  to  mass  flux  of  the  air  were 
investigated  over  the  range  of  O.O98  to  6.35  gm/cm2-sec  for 
combustor  diameters  of  1.?  to  7.62  centimeters. 

The  primary  tests  for  combustion  efficiency  and  temp- 
erature were  accomplished  using  1.9-centimeter  diameter 
zlrconia  tubss.  The  other  combustors  were  used  to  evaluate 
the  ease  of  initiating  steady  state  combustion,  the  effects 
of  combustor  heat  loss,  the  changes  in  th*  wall  deposits 
with  air  mass  flux  and  the  effect  of  combustor  length  on 
combustion  efficiency.  The  listing  of  the  combustors  and 


test  parameters  is  as  follows: 


Combustor  Material 


(/°u) 


'cm2-sec! 


Length 

(cm) 


Thin-Walled 
Stainless  Steel 


Graphite 


Zlrconia 


2.54  0.488-3.42  0.  -40 

5.08  0.195-0.488  20.  -61 

7.6 2 0.098-0.195  20.  -76 

2.54  0.977  15. 2-61 

1.9  1.71  -6.35  . 22.8-56 

2.54  0.488-3.42  30.5-76 
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Combustor  Materials 

The  tubular  steel  chambers  were  fabricated  from  24? 
stainless  steel.  The  graphite  burner  was  made  from 
stabilized  graphite  by  the  University  of  Dayton  under  a 

i 

contract  for  the  AP  Macerials  Laboratory.  It  consisted  of 

j 

7.62-centimeter  CD  graphite  cylinders  with  a 2.54-centimeter  j 

inner  core.  The  graphite  cylinders  were  Incased  in  a 7.62- 
centimeter  pipe  and  held  in  place  by  end  flanges;  the  exit 
end  flange  was  spring  loaded  to  allow  for  thermal  expansion. 

The  zirconia  chambers  were  formed  with  water  base 
zirconia  cement.  Glass  wool  fibers  were  mixed  into  the 
cement  to  provide  a matrix  binder.  Waxed  paper  was  wrapped 
around  a length  of  metal  tube  to  form  the  Inside  diameter  of 
the  combustor  and  the  plaster  mixture  was  trowelled  around  • 

the  metal  tube  in  layers  of  approximately  0.6  to  1.0- 
centimeter  thicknesses  until  a five-centimeter  thick  chamber 
wall  was  built  up.  Each  layer  was  allowed  to  cure  at  room 
temperature  before  the  next  layer  was  applied.  After  the  1 

plaster  had  cured,  the  metal  tube  was  removed  and  a steel 

t 

tube  was  inserted  into  the  combustor  entrance  to  enable  the  1 

combustor  to  be  connected  to  the  upstream  tubing.  The  waxed  i 

I 

5 

paper  burnt  away  immediately  upon  applying  the  initial  pre-  j 

heat  leaving  only  an  insulated  combustor.  j 

j 

Fuel  System 

1 

A fuel  system  was  developed  to  provide  a method  for 
obtaining  a steady  uniform  flow  of  fuel  particles  and  measur- 
ing the  fuel  flow  rate.  This  system  consisted  of  amorphous 

\ j 
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boron  pressed  Into  cylindrical  pellets  and  fed  into  a 
particle  mill  at  a known  rate. 

Peed  System.  The  feed  system  was  composed  of  a tubular 
fuel  chamber  and  a worm  gear  driven  piston.  The  worm  gear 
was  driven  by  a one  hp  variable  speed  electric  motor.  The 
speeds  of  the  piston  throughout  its  range  of  operation  were 
checked  periodically  during  the  experimental  investigation. 

The  maximum  deviation  of  the  speeds  at  any  motor  potentiom- 
eter setting  was  5#  with  the  root  mean  square  deviation  less 
than  I#.  The  rate  of  fuel  flow  was  calculated  from  the 
speed  of  the  piston  drive  multiplied  by  the  density  per  unit 
length  of  the  fuel  pellets. 

Particle  Mill.  The  particle  mill  consisted  of  an 
aluminum  chamber  and  a row  of  circular  wire  brushes,  as  shown 
in  Fig  6.  Three  Osborne  ringlock  No.  2679  brushes  with  a 
5.3-centimeter  diameter,  ll-M diameter  steel  wire  and  a 1.525- 
centimeter  arbor  hole  were  used.  The  brushes  were  rotated  at 
a constant  speed  of  6000  RPM.  The  aluminum  chamber  was 
constructed  so  that  there  were  a minimum  of  void  spaces  for 
dust  accumulation.  A chevron  seal  around  the  drive  shaft 
precluded  the  loss  of  fuel  particles  and  air.  To  determine 
the  effectiveness  of  the  particle  mill  returning  the 
amorphous  boron  back  to  its  original  particle  size,  a sample 
of  the  amorphous  boron  as  it  was  received  and  a sample  of 
repowdered  particles  were  analyzed  for  size  distribution. 

A Stokes  flotation  analysis  of  these  samples  was  performed 
by  the  Chemistry  Laboratory  of  the  AP  Aerospace  Research 
Laboratories.  The  results  from  their  analysis  is  suown  in 
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SIDE  VIEW 


A Fuel  Pellet  Feed  Chamber 
B Air  Flow  into  the 
•Particle  Mill  Chamber 
C Particle  Laden  Air  Flow 
to  Combustion  Chamber 
D Wire  Brushes 
E Arbor  Drive 


Figure  6.  Boron  Particle  Mill 
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Fig  7 as  the  aecumulat  fraction  of  the  sample  mass  with 
particle  diameters  less  than  that  indicated  on  the  abscissa. 
The  repowdered  boron  had  a smaller  size  distribution  than 
dl^  the  original  amorphous  boron.  Assuming  spherical 
particles,  the  number  of  fuel  particles  within  a diameter 
class  was  determined  from  the  mass  of  fuel  in  each  class. 

The  average  particle  diameter  was  then  calculated  from  this 
Information.  The  calculated  average  particle  size  for  the 
repowdered  boron  was  1.25  At  while  the  original  amorphous 
powder  had.  an  average  particle  size  of  1.75/a. 

Fuel  Preparation.  The  fuel  pellets  were  fabricated  by 
mixing  tap  water  with  Trona  amorphous  boron  powder  until  u 
thick,  smooth  paste  was  obtained.  The  paste  was  placed  in 
a teflon  lined  mold  and  pressed  slowly  on  a metallurgical 
press  to  6000  psi  to  remove  the  excess  water  and  boron 
trioxide.  The  pressure  was  applied  slowly  to  equalize  the 
pressure  within  the  mold  to  obtain  a homogeneous  density. 

The  pellets  were  extruded  from  the  mold  and  allowed  to  dry 
at  the  room  atmosphere.  After  they  were  thoroughly  dried, 
they  were  weighed  on  a chemical  balance  scales  and  measured 
to  obtain  their  linear  density.  The  residual  moisture 
content  of  a sample  pellet  was  checked  by  the  Fuels  Branch, 
AF  Aero  Propulsion  Laboratory,  Wright-Patterson  AFB,  Ohio, 
and  was  found  to  be  that  of  the  surrounding  atmosphere. 

Igniter  Flame.  An  igniter  tube,  constructed  of  0.64-cm 
stainless  steel  tubing,  was  located  at  the  aft  end  of  the 
combustor  aligned  parallel  with  the  combustor  center  line 

as  shown  in  Fig  5.  > Hydrogen,  passed  through  this  tube,  was 

»* 
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used  to  preheat  the  combustor  walls  and  to  Ignite  the  boron 
particles.  The  test  procedure  was  to  Ignite  the  hydrogen 
flame  at  a reduced  air  flow  then  stabilize  the  air  mass  flow 
to  the  desired  test  condition.  When  the  air  flow  was 
stabilized,  the  boron  particle  flow  was  introduced.  Once 
the  boron  A'lame  was  established,  the  hydrogen  was  turned  off. 
Ins  trumenta tlon 

Air  Flow.  A GL-8  Cox  turbine  type  flow  meter  was  used 
to  measure  the  air  flow  rate.  The  flow  meter  was  calibrated 
at  the  gas -.flow  facilities  at  Newark  APS,  Ohio  and  Tinker 
AFB,  Oklahoma  against  a known  air  volumetric  flow  at  20,  30 
and  40  psl.  The  meter  output  was  converted  to  an  analog 
signal  and  displayed  on  a Brown  recorder.  The  air  pressure 
in  the  constant  pressure  chamber  was  adjusted  to  the  desired 
test  condition  and  monitored  on  a calibrated  Heise  pressure 
gauge. 

Fuel  Flow.  A variable  speed,  one  hp,  electric  motor 
was  used  to  drive  the  fuel  piston.  The  motor  RPM  was  con- 
trolled by  a calibrated  potentiometer.  The  piston  drive 
speed  and  potentiometer  settings  were  correlated  by  a stop 
watch.  • 

Temperature  Measurement.  The  luminous  flame  emissions 
were  scanned  using  a Wamer-Swasey  501  Rapid  Scan  Spectrom- 
eter, both  in  scanning  mode  and  as  a monochrometer.  In 
the  scanning  mode,  the  spectrum  from  0.4  to  0.65  fk  was 
monitored.  The  boron  oxide  emission  peak  at  0.518  jk  was 
monitored  monochromatlcally.  The  spectrometer  signal  was 
enhanced  through  photomultipliers  and  displayed  on  an 
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oscilloscope.  Absorption  d;  ta  of  the  flame  were  obtained 
by  viewing  a calibrated  tungsten  lamp  through  • le  flame. 
Absorption-emission  information  was  obtained  by  chopping 
the  light  from  the  tungsten  lamp.  The  data  were  recorded 
on  Type  52  Polaroid  film.  The  method  of  calculating  the 
flame  temperatures  from  this  data  is  given  in  Appendix  G. 

Exhaust  Gases  Measurement.  Exhaust  gas  samples  were 
collected  through  a vacuum  operated  0.64-centimeter 
stainless  steel  probe.  The  samples  were  trapped  in  glass 
sample  bottles- with  high  vacuum  valves.  Oxygen/nitrogen 
and  nitrogen/argon  ratios  were  determined  by  the  High 
Energy  Chemistry  Branch  at  the  AP  Aerospace  Research 
Laboratories  using  mass  spectroscopy. 
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IV.  RESULTS  AND  DISCUSSION  j 


Numerous  test  runs  were  performed;  some  from  which  data 
were  collected  and  others  were  made  to  observe  the  appearance 
of  the  flame,  the  combustor  wall  deposits  or  the  effects  of 
the  combustion  process  on  the  combustor  materials.  The  test 
runs  were  reproducible  and  the  technique  for  establishing 
the  self -sustained  flame  became  apparent.  Self-sustained 
flames  of  particulate  boron  and  air  could  not  be  established 
for  any  test  condition  without  boron  reacting  on  the  combus- 
tor walls..  The  difficulty  of  establishing  steady  flames 
Increased  'as  the  fuel/air  ratio  decreased  and  at  equivalence 
ratios  less  than  0.5*  a steady  self -sustained  flame  could 
not  be  established.  Although  the  equivalence  ratio  at  0.7 
could  be  established  as  a steady  self -sustained  flame  with 
complete  burning  of  the  boron,  the  emissivity  of  a 
2.54- centimeter  diameter  flame  was  too  weak  and  lrratlc  to 
measure  the  temperature.  As  the  fuel/air  ratio  was  increased, 
the  flames  were  Increasingly  easier  to  establish  with  the 
partially  burned  boron  particles  of  fuel  rich  flames  burning 
to  completion  in  the  exhaust  plume.  But  above  an  equivalence 
ratio  of  three,  the  white  exhaust  cloud,  characteristic  of 
complete  combustion,  became  dark  brown  to  black  which 
indicated  many  of  the  particles  did  not  ignite  and  unburned 
boron  was  present  in  the  exhaust  plume. 

The  combustion  process  investigated  here  was  for  one 
fuel  particle  size  distribution  (with  an  average  diameter 
near  one /x ) and  no  Information  was  derived  for  the  influence 
of  different  particle  diameters.  However,  the  solids 
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loading  In  the  gas-particle  mixture  was  varied  to  determine 
the  effect  of  the  fuel/air  equivalence  ratio  on  the  combus- 
tion characteristics.  The  data  for  the  oxygen  content  and 
temperature  of  the  exhaust  gases  at  the  combustor  tube  exit 
are  tabulated  in  Tables  IV,  V and  VI  of  Appendix  J for 
equivalence  ratios  1,  1.5  and  2.  These  data  were  obtained 
from  1.9-centimeter  and  2.5^-centimeter  diameter  combustors 
with  lengths  varying  from  15  to  76  centimeters  and  are 
portrayed  in  Figs  8 through  27.  They  indicated  the  results 
wore  reproducible  and  a steady,  self -sustained,  boron  powder 
combustion  flame  could  be  established. 

The  discussion  is  arranged  in  three  parts:  ia)  the 

experimentally  measured  characteristics  are  presented  to 
show  that  this  combustor  design  technique  for  boron/alr 
combustion  resulted  in  rapid  complete  reaction  with  a high 
thermal  energy  release,  (b)  these  characteristics  are 
compared  to  the  analytical  combustion  model  of  this  study  to 
illustrate  the  relative  importance  of  the  heat  producing  and 
transfer  mechanisms  and  the  empirically  determined  character- 
istic values,  and  (c)  discussion  of  some  characteristics  of 


the  combustion  process. 

Fuel  Pellets 

The  fuel  pellets  obtained  from  the  fuel  fabrication 
process  were  1.75  centimeters  in  diameter  with  density  near 
0.9  gm/cc.  This  procedure  resulted  in  uniform  density 
pellets.  Approximately  10#  of  the  pellets  were  weighed  and 
measured,  these  had  an  rms  deviation  of  less  than  2#  from  the 
average  among  the  pellets.  Most  of  this  variation  was 
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caused  by  the  compressibility  of  the  press  mold  teflon  liner 
which  allowed  the  pellet  diameter  to  vary  slightly. 

Gas-Particle  Mixture 

The  uniformity  of  the  gas-particle  flow  was  checked  by 
measuring  the  attenuation  of  the  intensity  of  the  calibrated 
tungsten  lamp.  The  tungsten  lamp  was  viewed  through  the  gas- 
particle  flow  by  the  Warner-Swasey  Rapid  Scan  Spectrometer 
as  shown  in  Pig  5«  Oscilloscope  traces  of  the  wave  length 
spectrum  from  0.4  to  0.65/*  were  photographed  with  ^-second 
exposure.  This  was  accomplished  for  particle  loadings 
corresponding  to  fuel/air  equivalence  ratios  of  0.1  to  1.0 
with  no  broadening  of  the  oscilloscope  trace.  This 
indicated  that  the  gas-particle  flow  was  steady  for  all 
particle  loadings  checked.  The  average  diameter  of  the 
particles  in  the  gas-particle  mixture  was  calculated  using 
the  absorption  coefficient  developed  in  Eq  9»  it  was  deter- 
mined to  be  near  | micron  for  the  proportionality  constant 
Z equal  to  one.  This  compares  well  with  the  1.25 /*  average 
diameter  calculated  by  the  AP  Aerospace  Research  Laboratory 
from  a Stokes  flotation  analysis. 

Combustion  Efficiency 

The  combustion  efficiency  for  the  gas-particle  mixture 
was  evaluated  by  measuring  the  change  in  the  oxygen  content 
of  the  carrier  gas.  This  was  accomplished  by  obtaining  the 
ratio  of  nitrogen  to  oxygen  from  macs  spectrographic  data  of 
the  cooled  exhaust  gases  and  comparing  it  to  the  nitrogen/ 
oxygen  ratio  of  air.  This  method  assumed  there  were  no 
nitrogen  reactions  of  any  consequence  and  the  primary 
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oxygen  reactions  were  with  the  boron.  To  assure  that  this 
assumption  was  true,  nitrogen/argon,  NO/argon  and  NOg/argon 
ratios  of  the  mass  spectrographic  data  were  made.  There  was 
no  change  in  the  nitrogen/argon  ratio  as  shown  in  Table  VII 
of  Appendix  J,  where  the  values  obtained  for  this  ratio  were 
within  the  natural  scatter  of  the  data.  The  data  for  NO  and 
NOg  were  too  small  to  be  noted.  Thus,  assuming  no  nitrogen 
reactions,  the  exit  oxygen  depletion  parameter  was 


The  validity  of  any  gas  sampling  data  is  always  in 
question.  The  major  reason  for  this  is  the  quenching  of 
the  sample.  For  solid  particles  reacting  in  a gas,  the 
rate  of  quenching  would  be  slower  than  for  a gaseous 
reaction.  The  sampling  procedure  for  a particle  system  was 
further  complicated  since  the  particles  were  either  molten 
or  had  molten  oxides  or.  their  surfaces.  Th^se  particles 
adhered  to  the  gas  probe  causing  restrictions  of  the  flow 
into  the  probe  or  deflection  of  the  exhaust  plume  away  from 
the  probe.  The  Internal  restrictions  caused  by  the  solids 
collecting  within  the  probe  were  reduced  by  adding  a helium 
purge  gas  system  to  the  sampling  probe.  These  solids  were 
removed  by  purging  the  probe  after  the  sample  was  taken. 

The  sampling  time  to  ensure  a represer.  tive  sample  of  the 
exhaust  gases  was  reduced  by  purging  the  probe  just  prior  to 
taking  a sample.  R&pid  quenching  of  the  boron  reaction  was 
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believed  to  take  place  because  boron  reacts  vigorously  only 
above  1800  to  2000°K  with  almost  no  reaction  below  1000°K. 
Thus,  the  oxygen  content  in  the  sample  was  assumed  to  be 
approximately  the  same  as  that  of  the  exhaust  gas  from  which 
it  was  obtained. 

Correlation  of  the  Exhaust  Gas  Data 

The  exhaust  gas  data  were  collected  from  different 
length  combustors  and  at  various  mass  fluxes  of  the  carrier 
air.  The  combustor  length  was  measured  from  the  pilot 
burner  tip. to.  the  combustor  exit  and  the  air  mass  flux  was 
calculated ' from  the  total  air  flow  into  the  system.  The 
oxygen  depletion  parameter  values  y were  plotted  against 
the  combustor  residence  time  parameter  4T  , defined  in  Eq  31, 
to  determine  the  trends  for  the  combustion  efficiency,  as 
shown  in  Figs  8 and  9.  The  data  at  4>  * 1.5  was  so  similar 
to  that  at  = 2.0,  it  was  plotted  with  the  data  for  2.0. 

At  low  air  mass  fluxes,  the  oxygen  depletion  data 
collected  were  shifted  toward  the  longer  residence  times. 

This  gave  the  indication  of  lower  reaction  rates  for  the  low 
gas  flows.  However,  this  shift  of  the  data  was  not  caused 
by  a lower  reaction  rate  but  was  due  to  excessive  deposits 
of  boron  and  its  oxides  in  the  upstream  end  of  the  combustor. 
These  deposits  filled  in  the  volume  at  the  entrance  of  the 
combustion  chamber  and  created  a shortened  combustor.  This 
tendency  of  the  particles  to  agglomerate  at  the  entrance  of 
the  combustor  decreased  as  the  air  mass  flux  increased  until 
above  approximately  1.5  gm/cm2-sec  the  wall  deposits  were 
uniform.  Thus,  the  majority  of  the  test  runs  recorded  were 
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near  3 gm/cm2-sec  to  insure  the  uniformity  of  the  wall 
deposits. 

Data  were  collected  only  for  tests  where  a stable 
luminous  flame  was  visible  at  the  exit  plane.  For  this 
condition,  the  minimum  oxygen  depletion  value  was  0.1 
at  - 1*0  and  approximately  0.4  for  equivalence  ratios  of 
1,5  and  2.0.  This  corresponded  to  a value  of  the  residence 
time  parameter  of  8.6  as  listed  in  the  tables  of  data  for  a 
1.9-centimeter  diameter  chamber,  22.8  centimeters  in  length. 

With  increased  residence  time,  90  to  100#  of  the  oxygen  was 
consumed  for  residence  time  parameter  values  greater  than 
twenty  milliseconds.  Even  though  the  higher  values  of 
oxygen  depletion  were  obtained  more  rapidly  at  an  equiva- 
lence ratio  of  2.0,  the  oxygen  depletion  rates  from  the  onset 
of  the  luminous  flame  to  complete  consumption  of  the  oxygen 
were  very  similar  for  both  equivalence  ratios  of  1.0  and 
2.0.  The  greatest  deviation  in  the  oxygen  depletion  was  at  \ 

B,  of  13.4  for  an  equivalence  of  one.  This  was  at  the 
steepest  slope  for  the  reaction  and  a slight  change  in  the 
flame  front  would  account  for  most  of  the  deviation. 

Correlation  of  the  oxygen  depletion  data  with  the  real 
time  combustor  residence  time  indicated  that  the  stable 
luminous  flame  was  established  in  approximately  two 
milliseconds  and  the  reaction  was  complete  in  four 
milliseconds.  This  was  evaluated  from  data  at  an  equiva- 
lence ratio  of  2.0  in  Fig  10;  the  times  would  not  be 
significantly  different  for  equivalence  ratios  of  1.0  r.  t 

1.5*  The  combustor  residence  time  was  calculated  from  the 
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residence  time  parameter  by  integrating  Eq  31  with  the  tem- 
pera ture  values  obtained  from  Pig  21  using  straight  line  seg- 
ments connecting  the  average  values  of  the  temperature  data. 
Temperature  Measurements 

Because  the  flame  temperatures  encountered  were  much 
too  high  for  thermocouple  measurements  and  due  to  the  solids 
in  the  flame  zone,  a noninterference  measurement  was  used. 

The  boron  flames  are  highly  luminous,  therefore,  an  absorp- 
tion-emission measurement  system,  described  by  Beyermann 
(fief  2?),.  was  adopted.  This  system  utilized  Planck*s  law  of 
radiation  and  is  theoretically  applicable  only  to  gaseous 
media. 

The  theoretical  methods  for  applying  the  absorption- 
emission  technique  to  temperature  calculation  are  shown  in 
Appendix  G.  The  flame  temperature  at  the  combustor  exit 
was  calculated  using  Eq  G-4.  This  technique  required  the 
monochromatic  radiant  intensity  of  a known  temperature 
source,  the  emitted  intensity  of  the  test  flame  and  a 
measure  of  the  attenuation  of  the  known  source  through  the 
test  flame.  The  monochromatic  absorptivity  of  the  flame  is 
proportional  to  the  difference  in  the  alternating  absorption- 
emission  intensity  signals  while  the  flame  intensity  is 
proportional  to  the  lower  (emission)  intensity.  Fig  29* 

Typical  absorption-emission  data  for  equivalence  ratios 
of  1.0  and  2.0  at  0.518  Jx.  are  shown  in  Pigs  11,  12  and  13, 
where  the  emitted  intensity  of  the  flame  is  the  lower  signal 
and  the  flame  plus  the  intensity  of  a controlled  tungsten 
lamp  viewed  through  the  flame  is  the  upper  signal.  The 
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Figure  11.  Spectrometric  Data  of  the  Flame  Emissions 
from  1.9-cm  ID  Zirconia  Combustor  at 
* = 0.518/4  for  <f>  = 2 at  f = 8.62  ms 
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Figure  12.  Spectrometrlc  Data  of  the  Flame  Emissions 
from  1.9-cm  ID  Zlrconla  Combustor  at 
>.  = 0.518/*  for  <f>  = 2 at  £ = 21.1  ms 


Figure  13.  Spectrometric  Data  of  the  Flame  Emissions, 
from  1,9-cm  ID  Zirconia  Combustor  at 
*=  0.518/4.  for  j>  = 1 at  § = 34.1  ms 
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chopped  signal  from  the  tungsten  lamp  provided  the  known 
temperature  calibration  scale.  As  can  be  seen  from  Pigs  11 
through  13»  the  signal-to-noise  ratio  was  quite  high 
especially  for  equivalence  ratios  near  unity.  Equivalence 
ratios  less  than  one  produced  signals  that  were  too  weak  and 
irratlc  to  deduce  a temperature.  The  noise  level  of  the 
flame  signal  also  increased  with  increasing  temperature. 
There  was  very  little  broadening  of  the  signal  between  the 
emission  and  absorption  scans;  this  indicated  that  the 
monochromatic,  absorptivity  of  the  flame  was  constant. 

To  calculate  a temperature  from  this  data,  there  was 
uncertainty  as  to  what  portion  of  the  signal  should  be  used 
to  evaluate  the  flame  eraissivity  and  intensity.  Flame 
temperatures  were  calculated  from  the  spectrometrlc  data 
using  various  positions  within  the  broadened  signals  to 
evaluate  the  change  in  calculated  temperature.  The 
uncertainty  in  the  calculated  temperature  was  proportional 
to  the  flame  temperature  and  inversely  proportional  to  the 
equivalence  ratio.  The  maximum  variation  in  a calculated 
flame  temperature  was  50°K  at  approximately  2200°K  for  an 
equivalence  ratio  of  two  and  over  1000°K  for  a flame  at  an 
equivalence  of  one  near  3&00°K. 

During  the  test  runs  at  high  combustion  efficiency  for 
an  equivalence  ratio  of  one,  the  zirconia  chamber  would 
become  so  molten  that  sections  of  the  chamber  wall  were 
blown  from  the  combustor.  It  was  only  at  an  equivalence  of 
one  that  this  condition  occurred.  After  these  tests,  the 
combustor  wall  near  the  exit  had  a layer  of  light  gray 
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material  on  yellow  zirconla  ceramic,  as  represented  by  the 
sketch  in  Fig  14.  The  gray  material  could  only  be  melted 
by  an  oxyacetylene  torch  adjusted  to  its  hottest  flam? 

(near  J000°K)  and  impinging  directly  on  a specimen  of  the 
substance.  These  materials  were  analyzed  by  the  Air  Force 
Materials  Laboratory.  The  bright  yellow  ceramic  was  identi- 
fied as  cubic  zirconla  and  the  light  gray  material  as  pure 
monoclinic  zirconla. 

The  melting  point  of  monoclinic  zirconla  was  reported 
in  Table  XI  of  Appendix  J to  be  2900°K  and  a flame  tempera- 
ture of  2940°K  was  calculated  for  the  first  maximum 
stoichiometric  test  from  absorption-emission  data  using  the 
maximum  signal  peaks.  With  this  information  as  a calibration 
point,  all  the  temperature  Information  was  calculated  using 
the  maximum  signals. 

The  cause  for  the  fluctuations  in  the  flame  signal  in 
Figs  11  through  13  is  not  known.  It  could  arise  from 
several  sources,  the  most  probable  are:  blocking  of  the 

flame  emissions  by  the  condensing  oxides  on  the  periphery  of 
the  flame,  scattering  of  the  signal  by  solids  in  the  flame, 
more  dense  oxides  surrounding  the  particles,  or  natural 
fluctuations  from  one  emission  peak  to  another..  The 
frequency  of  the  fluctuation  increased  and  the  amplitude 
decreased  with  increasing  equivalence  ratios;  thus,  the 
unburned  boron  or  oxides  surrounding  the  particles  could 
cause  the  fluctuations.  However,  the  greatest  amplitude  of 
the  fluctuations  occurred  at  high  combustion  efficiency  at 
an  equivalence  of  one  where  all  the  fuel  particles  should 
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have  burned  up  and  all  constituents  of  the  exhaust  plume 
should  have  become  gaseous.  The  boron  oxide  absorption- 
emission  spectra  for  the  observed  flame  signals  is  described 
in  the  literature  as  "fluctuation  bands"  without  explanation 
of  their  definition  of  fluctuation  (Refs  23»  24  and  25)* 

Figure  15  shows  a typical  spectra  of  the  boron/air 
flame  at  an  equivalence  of  one  and  the  emission  lines  are 
identical  to  those  described  in  the  literature.  This 
spectra  was  taken  in  the  rapid  scan  mode  of  the  Warner-Swasey 
Spectrometer  at  two  milliseconds  per  scan;  the  exposure  time 
was  twenty  milliseconds.  The  lower  signals  were  the  flame 
emissions  and  the  upper  signals  were  of  the  flame  and  the 
reference  tungsten  lamp.  The  emission  fluctuations  were 
clearly  illustrated  here  and  the  steady  nature  of  the  fluc- 
tuations was  shown  in  Figs  11*  12  and  13  where  the  0.518 JJL 
line  was  scanned  at  one  centimeter  per  second.  The  maximum 
emission  signal  was  produced  by  the  line  at  approximately 
0. 518yd  . As  can  be  seen,  this  maximum  emission  signal  was 
much  greater  than  the  continuum  emissions;  thus,  the  gaseous 
emissions  were  the  dominant  portion  of  the  signal.  The 
oxides  condensing  on  the  periphery  of  the  plume  would  be 
much  cooler  than  the  core  of  the  exhaust  flow.  If  this  was 
the  primary  cause  of  the  fluctuations,  then  the  best  explana- 
tion of  the  flame  signal  would  be  that  the  hot  core  of  the 
flame  was'  seen  as  the  maximum  portion  of  the  signal  and  would 
best  describe  the  exhaust  temperature. 

The  change  in  the  exhaust  temperature  with  increasing 
combustor  residence  time  was  correlated  to  the  combustion 
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Figure  15.  Typical  Oscilloscope  Trace  at  Exha  ;st  Plane 
Showing  a Characteristic  Boron  Flame  Spectra 
for  Maximum  Flame  Temperature  at  6 - 1 
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model  from  this  study,  Eqs  33  through  38,  for  an  equivalence 
of  one  in  Fig  1 6 and  for  an  equivalence  of  two  in  Fig  17 
with  the  residence  time  parameter  £ . The  calculated 
temperature  data  at  the  different  values  of  the  residence 
time  parameter  were  taken  from  zirconia  combustors  with  the 
same  lengths  and  air  mass  fluxes  and  thus  represented  the 
spread  of  the  data.  The  minimum  temper  rnrcs  measured  were 
from  2100<>K  to  2200°K.  These  were  not  the  minimum  tempera- 
tures for  the  luminous  flames  but  were  the  minimum  exhaust 
temperatures  .for  a steady,  self-sustained  flame  to  become 
stabilized’ in  the  tubular  combus to:--"! , The  fluctuation 
characteristics  of  the  .lame  signal  in  Fig  11  indicated  that 
that  particular  flame  was  not  stable.  The  temperature 
continued  to  Increase  as  the  residence  time  parameter  was 
increased  with  the  maximum  temperature  approaching  the 
adiabatic  flame  temperature  for  as  shown  in  Fig  18. 

Further  .trends  in  the  temperature  increase  for  <j>  - 2 in 
Fig  17  indicated  that  the  tempera  cure  rose  rapidly  to 
asymptotically  approach  the  B2O3  adiabatic  flame  temperature « 
With  further  Inert- a;  e in  combustor  length  and  well  above 
combustor  lengths  for  complete  oxygen  depletion,  the  temper- 
ature began  a second  increase  above  the  maximum  temperature 
for  BgO y This  econd  rise  indicated  a possible  boron 
nitride  reaction;  however,  due  to  the  large  source  of 
nitrogen  in  the  air,  this  was  not  substantiated  by  a change 
in  the  Ng/argon  ratio.  This  trend  was  duplicated  for 
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<f>-  1.5*  in  Table  V of  Appendix  J,  but  was  not  present  in 
the  stoichiometric  fuel/air  ratio. 

Maximum  Flame  Temperatures 

The  flame  temperatures  for  test  runs  at  high  combustion 
efficiency  or  runs  where  the  combustor  residence  time  was 
more  than  sufficient  to  yield  95  to  100^  oxygen  depletion 
were  plotted  versus  equivalence  ratio  in  Fig  18.  These 
temperature  data  were  collected  from  the  zirconia  combustors 
where  the  combustion  process  was  near  adiabatic.  The 
experimentally-  determined  temperatures  were  compared  to 
theoretical  predictions  for  adiabatic  temperatures  obtained, 
at  equilibrium  conditions  and  for  boron  trio.cide  as  the  only 
product.  Some  of  the  temperatures  deduced  from  the  data 
at  <j>  = 1.0  and  1.5  fell  between  the  two  curves  indicating 
that  the  chemical  heat  release  was  greater  than  that 
predicted  by  the  equilibrium  process.  Since  the  maximum 
measured  temperatures  at  these  equivalence  ratios  were 
greater  than  the  adiabatic  flame  temperature  for  equilibrium 
but  less  than  the  adiabatic  flame  temperature  for  BgO^,  the 
preferred  reaction  product  appeared  to  be  B20-j. 

Combustion  Model 

The  next  sections  of  this  discussion  examine  the 
phenomena  observed  during  the  experimental  investigation 
that  verify  the  theoretical  model  formulated  in  chapter  two 
and  evaluate  the  parameters  identified  in  the  theoretical 
analysis.  The  phenomena  were  the  role  of  the  heat  flux  from 
the  flame  in  the  heat  up  of  the  gas-particle  mixture;  the 
importance  of  the  combustor  wall  reactions  and  their 
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attendant  radial  heat  flux;  and  the  contribution  to  the  flame 
absorption  coefficient  due  to  the  combustion  products.  The 
parameters,  listed  on  page  39,  chapter  two,  were  evaluated 
by  correlating  the  experimental  results  with  the  theoretical 
analysis.  The  results  for  the  theoretical  model  were 
arrived  at  by  solving  Eqs  33  through  38  using  a Runga-Kutta 
technique.  The  associated  Fortran  Program  is  listed  in 
Appendix  H. 

Temperature  and  oxygen  depletion  parameter  data  for 
<j>  of  1,0,.  1,5  and  2.0  were  correlated  in  Fig  19;  the 
adiabatic  ‘tempera ture  for  B20^,  as  calculated  from  the  JANAF  _ 
tables,  was  included  as  a comparison.  At  the  low  oxygen 
depletion  values,  the  measured  temperatures  were  greater 
than  those  which  could  be  obtained  by  the  adiabatic  process 
indicating  that  external  energy  was  supplied.  At  the  higher 
oxygen  depletion  values,  the  temperature  deviation  was  less 
pronounced  with  the  higher  equivalence  ratios  coinciding 
with  the  adiabatic  temperature  trend;  the  measured 
stoichiometric  temperatures  continued  to  approach  the 
adiabatic  values  with  the  maximum  measured  temperatures  less 
than  the  maximum  adiabatic  temperature. 

These  temperatures  indicated  that  energy  to  heat  the 
gas-particle  mixture  above  the  adiabatic  values  was  supplied 
by  the  heat  transferred  from  the  flame.  The  mass  flux  of 
air  necessary  for  the  flame  to  supply  all  of  the  energy 
required  to  heat  the  oncoming  particles  to  ignition  was 
calculated  using  the  relation  developed  by  Essenhigh  and 
Csaba  (Ref  19) 
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When  Tp  = 2500°K,  T^  = 1850°K  and  £p  = 1. 

This  value  of  the  mass  flux  of  air  would  be  the  maximum 
possible  since  it  assumed  the  flame  to  be  a black  body. 

This  value  would  be  increased  by  a factor  of  two  if  the 
adiabatic  flame  temperature  at  ^ of  1.0  was  used;  however, 
the  mass  flux  of  air  would  still  be  a factor  of  ten  less 
than  that  used  in  this  investigation.  Therefore,  it  was 
concluded  that  heat  addition  from  boron  reacting  on  the 
combustor  wall  was  required  to  stabilize  the  flame. 

Radiation  Stabilized  Flame 

Tests  were  made  to  determine  whether  the  radiation  from 
the  luminous  flame  could  support  the  gas-particle  mixture 
heat  up. to  ignition  temperature  without  the  wall  reactions. 

A boron  flame  with  a hydrogen  pilot  flame  was  established 
in  a 7.62-centimeter  steel  chamber  at  an  air  mass  flux  of 
0.1  gm/cm2-sec  with  a fuel  particle  loading  for  an  equiv- 
alence ratio  of  two.  However,  a self-sustained  flame  could 
not  be  stabilized  without  the  presence  of  radiant  wall 
deposits.  The  combustor  wall  deposits  were  initiated  by 
preheating  the  walls.  Once  the  walls  were  coated  with  a 
radiating  layer  of  boron  and  its  oxides,  the  flame  emissivity 
increased  sufficiently  so  that  the  hydrogen  pilot  could  be 
turned  off.  The  boron  continued  to  be  heated  to  ignition. 

As  the  combustion  of  the  boron  continued  to  increase  in 
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intensity,  the  thickness  of  the  luminous  flame  increased, 
increasing  the  emissivlty  of  the  flame. 

With  increased  flame  luminosity,  the  flame  front  flashed 
back  through  the  oncoming  gas -particle  mixture  to  become 
stabilized  at  the  chamber  entrance.  This  demonstrated  that 
the  heat  flux  from  the  flame  could  cause  the  flame  to 
propagate  into  a boron  particle  cloud. 

Combustor  Wall  Reactions 


To  determine  the  characteristics  of  the  wall  deposits 
and  to  ascertain  that  they  were  primarily  burning  boron  and 
not  molten  oxides  heated  by  the  flame,  a series  of  tests 
were  made  without  a flame  being  established.  This  was 
accomplished  by  positioning  the  igniter  tube  near  the  exit 
of  a steel  chamber  at  an  air  mass  flux  of  l.?6  gm/cm2-sec 
and  an  equivalence  ratio  of  one.  A brilliant  green-white 
flame  was  established  in  what  would  be  the  hydrogen  flame 
but  it  would  extinguish  when  the  hydrogen  pilot  was  turned 
off.  The  igniter  tube  was  retracted  into  the  steel  tube 
until  a short  deposit  was  formed.  The  flame  would 
extinguish  as  soon  as  the  pilot  flame  was  removed  but  the 
deposit  layer  continued  to  glow  a bright  red  as  long  as  the 
powder  flow  was  sustained.  This  process  was  repeated  until 
a deposit  length  of  approximately  20  centimeters  was  formed 
The  deposit  would  continue  to  glow  under  a sustained  fuel 
flow  without  the  green-white  flame,  only  red  streaks  of 
partially  ignited  boron  were  observed  in  the  exhaust  dust. 
The  steel  tube  with  these  deposits  glowed  a bright  red  in 
color.  The  color  temperature  of  the  tube  was  observed 
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using  a pyrometer  with  the  measured  temperature  almost 
uniform  at  2000°F.  This  indicated  that  for  a uniform 
thickness  deposit,  the  heat  flux  from  the  reactions  inside 
the  tube  was  not  constant  along  the  tube..  By  increasing  the 
fuel  flow,  the  boron  cloud  could  be  ignited  and  burn  as  the 
characteristic  green-white  flame. 

Evaluation  of  the  Heat  Up  Zone  (Zone  I) 

The  values  of  oxygen  depletion  and  temperature  were 
measured  at  the  exit  of  different  lengths  of  combustors. 
However,  for  the  combustion  model,  it  was  assumed  that  the 
combustor  lengths  at  which  the  data  were  taken  corresponded 
to  the  various  locations  along  a single  combustor.  In  the 
combustion  model,  the  upstream  conditions  were  affected  by 
the  downstream  conditions.  The  same  downstream  conditions 
did  not  exist  for  the  shorter  combustors  where  the  reaction 
was  not  complete  but  were  replaced  by  a luminous  exhaust 
plume.  The  exhaust  plume  characteristics  as  seen  by  the 
oncoming  gas-particle  mixture  were  assumed  to  be  similar  to 
those  of  the  confined  flame  for  this  model. 

No  experimental  measurements  were  taken  in  the  heat  up 
region  of  the  gas-particle  mixture.  Nor  was  there  any 
attempt  to  measure  the  temperature  at  which  the  reaction  of 
the  particles  within  the  mixture  became  self-sustaining. 
However,  temperature  and  oxygen  depletion  parameter  data 
were  measured  near  the  minimum  residence  time  parameter 
necessary  to  sustain  a steady  luminous  flame.  The  para- 
metric values  of  the  average  flame  temperature  Tp,  the 
proportionality  factor  Z,  and  the  wall  reaction  mode] 
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(d'y//d§)w  were  evaluated  using  Eqs  33  through  37  to 
approach  these  measured  values.  These  parameters  were 
interrelated  where  an  assumption  of  a value  for  Z less  than 
unity  would  require  an  increase  in  the  values  for  Tp  or 
(d1^/d£  ) . Correlated  with  data  at  equivalence  ratios  of 

W 

1.0  and  2.0,  the  resulting  values  for  Z equal  to  one  were: 

Tp  equal  to  2900°K  for  <f>  of  1.0  and  2400°K  for  <f>  of  2.0 
and  (d|>/d£  )w  = 0.014  ms"1  or  g = 4.67  x 10“2  (sec-K)"1. 
These  results  were  compared  to  the  experimental  measurements 
in  Figs  20  .and  *21. 

The  absorption  coefficient  proportionality  factor  Z is 
a measure  of  the  particle  absorptivity  and  Z approximately 
equal  to  one  indicated  that  the  boron  particles  approach 
black  body  absorptivity.  This  value  of  Z was  the  same  as 
that  assumed  for  the  cold  gas-particle  mixture,  page  52. 

This  result  yields  an  absorption  coefficient  k near  unity  at 
room  temperature.  But,  the  absorption  coefficient  is  propor- 
tional to  the  air  density  and  decreases  as  the  mixture 
temperature  increases.  Therefore,  as  the  gas-particle 
mixture  approaches  the  self -sustained  temperature,  where  the 
major  radiant  energy  transfer  from  the  flame  is  applied,  the 
absorption  coefficient  is  much  less  than  unity  and  the 
optically  thin  limit,  assumed  for  this  model,  is  approached. 

The  equivalence  ratio  dependence  of  the  heating  due  to 
the  wall  reaction  and  the  radiant  energy  transfer  from  the 
flame  provided  an  explanation  for  the  difficulty  in  obtain- 
ing a self -sustained  flame  in  fuel  lean  or  extremely  fuel 
rich  mixtures.  The' energy  release  from  the  wall  reactions. 
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so  necessary  for  a self-sustained  flame  decreased  as  the 
equivalence  ratio  was  decreased.  Also,  the  adiabatic  flame 
temperature  shown  in  Fig  18  approaches  2200  to  2100°K, 
the  minimum  temperature  necessary  for  self-sustained  combus- 
tion, as  the  equivalence  ratio  is  decreased  below  0.5  or 
increased  above  3*0. 

Particle  Reaction  Model 

By  correlating  the  measured  oxygen  depletion  and  temper- 
ature data  with  Eq  37*  the  constants  A*  and  Eft/H  were 
evaluated...  The  value  for  the  activation  energy  is  quite 
insensitive  to  change  with  values  of  Ea/R  from  10^  to 
5 x 10^  °K  giving  similar  oxygen  depletion  rates.  These 
trends  and  the  dependence  of  on  Ea/R  are  shown  in  Figs 
22  and  23  for  A*  of  1423  sec”1.  The  correlations  were 
calculated  using  an  average  of  the  measured  temperatures 
with  Ea/R  ?f  2 x giving  the  best  correlation  to  the 
oxygen  data.  The  Arrhenius  pre-exponential  constant  A was 
then  evaluated  in  Appendix  I and  was  calculated  to  be 
9.04  x 10“3  moles  of  air/cm-sec. 

Combustion  Model  for  the  Flame  Zone  (Zone  II) 

Temperature  calculations  using  the  equations  for  the 
combustion  model  for  Zone  I,  Eqs  33  through  37 » were 
continued  after  the  cloud  particle  reactions  became  the 
dominant  source  of  energy.  This  resulted  in  calculated 
temperatures  well  above  the  measured  values.  Much  of  this 
was  due  to  the  form  of  the  radiant  heat  transfer  term  which 
continued  to  be  a heat  source  until  the  average  flams 
temperature  was  reached  and  did  not  provide  a mechanism  to 
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transfer  energy  to  the  oncoming  cloud.  The  rate  of  Increase 
in  measured  temperature  was  less  than  the  rate  of  increase 
in  adiabatic  temperature  at  the  higher  oxygen  depletion 
values,  as  seen  in  Pig  19»  indicating  an  energy  loss  from 
the  flame  zone  at  this  point.  The  boron  flame  is  very 
luminous  and  as  indicated  in  Pig  15 » the  luminosity  resulted 
from  continuum  radiation  as  well  as  from  radiating  gases. 

The  fuel  particles,  especially  for  the  fuel  rich  flames,  and 
the  reaction  products,  condensed  and  gaseous,  contributed  to 
the  flame,  luminosity. 

Equation  38,  along  with  Eqs  34,  35  and  37 t was  evaluat- 
ed for  the  flame  zone  to  determine  thr,  contribution  of  the 
reaction  products  to  the  flame  mass  absorption  coefficient 
Xg«  The  value  of  X g for  temperatures  less  than  2100°K 
was  considered  negligible.  At  2100°K,  X g was  a maximum 
indicating  the  transition  from  Zone  I to  Zone  II.  Its  value 
decreased  steadily  with  increasing  'Y'  to  become  negligible 
above  the  oxygen  depletion  value  of  0.8.  The  linearized 
relation  for  X use(*  for  flame  zone  calculation  in 

O 

Appendix  H to  correlate  the  experimental  data,  was  calculated 
to  be 


Xj  s °l  [ l+O.  (S-T/MOP) X /<?  (39) 

This  value  is  considerably  larger  than  the  particle  mass 
absorption  coefficient  X p used  the  correlation  which 
was 
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The  calculated  temperatures  and  oxygen  depletion  values 
were  compared  to  measured  data  in  Pigs  8,  9,  16,  17,  20  and 
21.  The  predicted  oxygen  depletion  values  at  equivalence 
ratios  of  1.0  and  2.0  and  temperature  trends  at  of  1.0 
agreed  with  the  measured  values  very  well.  The  calculated 
temperatures . compared  to  the  measured  values  in  Pigs  17  and 
21  for  an  equivalence  ratio  of  2.0  indicated  that  the  flame 
emlssivity  was  underestimated  at  2100°K  to  2200°K  and  over- 
estimated near  2300°K.  The  flame  emlssivity  was  considered 
negligible  above  ^ of  0.8  in  the  calculations  and  this 
resulted  in  an  overprediction  of  the  temperature  rise. 

The  sudden  increase  of  the  combustion  model  flame 
emlssivity  corresponded  to  the  rapid  rise  in  the  reaction 
rate  of  the  suspended  particles  and  the  decrease  in 
emlssivity  followed  the  reaction  rate  as  the  oxygen  was 
consumed.  The  measured  values  for  the  gaseous  absorption 
coefficient  at  = 0.518  yw.  followed  this  trend.  This 
result  is  shown  in  Pigs  24  and  25.  The  measured  values  of 
the  absorption  coefficient  at  0.518M  were  less  than  the 
values  calculated  using  Eq  39.  This  was  as  expected  since 
the  measured  values  were  at  one  wave  length  while  the 
calculated  values  were  for  the  entire  spectrum.  Both  the 
measured  and  calculated  absorption  coefficients  indicated 
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that  the  maximum  radiant  emissions  for  a boron  flame 
occurred  Just  above  the  minimum  temperature  necessary  for 
self -sustained  flames  (2100  to  2200°K). 

Gaseous  Emlsslvlty 

John  (Bef  24)  attributed  the  absorption-emission  lines 
from  0.4  to  0.6 to  B02  only  and  not  from  gaseous  B2O3. 
Thus,  the  value  of  the  absorption  coefficient  at  0.518 yu 
was  a measure  of  the  concentration  of  the  B02  species.  The 
gaseous  absorption  coefficient  at  0.518  ylt  followed  the 
particle  reaction  rate  closely  indicating  that  the  concentra 
. tlon  of  BO 2 was  proportional  to  the  reaction  rate  and  a.  „ 
possible  short-lived  intermediate  reaction  species.  A 
further  indication  of  this  was  the  correlation  of  the 
volumetric  absorption  coefficient  to  1/T  in  Pig  26.  The 
data  from  2200°K  to  3000°K  for  all  equivalence  ratios  were 
correlated  by  least  squares  for  a slope  equal  to  2 x 10^. 
This  slope  corresponded  to  the  empirical  activation  energy 
calculated  for  the  reaction. 

Combustor  Wall  Deposits 

The  growth  of  the  wall  deposits  and  the  stability  of 
the  self -sustained  flame  was  achieved  much  faster  when  the 
chamber  walls  were  preheated.  During  combustion,  the 
deposits  glowed  a brilliant  red-white  and  were  molten. 

After  the  deposits  oooled,  they  had  a very  hard,  brittle, 
porous  structure  with  a dark  brown  to  black  coloration. 
Samples  of  these  deposits  were  analyzed  by  the  Air  Force 
Materials  Laboratory  which  reported  the  major  species  to  be 
B203  with  traces  of'  BgO  and  BN. 
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The  thickness  of  the  deposits  appeared  to  be  dependent 
on  the  air  mass  flux  and  the  thermal  conductance  of  the 
combustor  materials.  The  thickness  of  the  deposit  layer 
was  greater  in  the  stainless  steel  chambers  and  the  least 
in  the  zirconia  combustors.  It  was  difficult  to  character- 
ize the  wall  deposits  in  the  graphite  chamber  as  they  did 
not  stabilize  but  continued  to  pour  out  of  the  burner  exit 
in  the  molten  state. 

For  the  zirconia  combustors  at  the  highest  air  mass 
flux  test  conditions,  the  entire  deposit  thickness  from  the 
upstream  end  to  the  exit  of  the  combustor  was  uniform  with 
a smooth  inside  surface.  As  the  air  mass  flux  was  decreased 
below  1.5  gm/om2-sec,  the  upstream  entry  section  deposit 
thickened  until  only  a very  small  entry  hole  remained. 

These  entry  deposits  in  the  zirconia  chambers  were  identical 
to  those  in  the  stainless  steel  tubes.  This  condition  was 
stable  and  the  tubes  did  not  completely  close  off  even  after 
numerous  repeated  runs.  Immediately  downstream  of  this 
small  entry  port,  the  deposits  thinned  rapidly  to  a minimum 
where  the  flame  was  fully  established.  This  thinner  deposit 
layer  occurred  as  the  self -sustained  flame  became  stabilized. 
At  the  beginning  of  each  repeated  run,  the  oxide  deposits 
would  become  molten  and  sluff  off  to  establish  the  steady 
combustor  conditions.  When  this  occurred  in  the  2,54- 
centimeter  stainless  steel  tubes,  it  resulted  in  the  failure 
of  the  tube  and  prevented  high  efficiency  combustion. 

The  dependence  on  the  thermal  conductance  was  also 
observed  in  the  ziroonia  chambers.  A layer  of  the  zirconia 
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plaster  under  the  boron  deposits  would  become  molten  and 
form  cubic  zlrconia  ceramic.  This  ceramic  layer  would 
thicken  during  repeated  runs.  The  boron  oxide  deposits 
were  approximately  0.15  centimeters  thick  after  an  Initial 
run  and  would  thicken  to  almost  0.6  centimeters  after  five 
to  ten  repeated  runs  for  a total  burner  life  of  ten  to 
fifteen  minutes  at  3 gm/cm2-sec  air  mass  flux.  This  change 
in  the  cross-sectional  area  of  the  combustor  due  to  the  wall 
deposits  was  not  apparent  in  the  data  collected.  Therefore, 
this  condition,  may  be  primarily  an  effect  of  the  cooled 
deposits  as  there  was  no  way  to  observe  the  deposits  during  . 
a combustion  run. 

Evaluation  of  the  Evaporation  Rate  Constant 

The  evaporation  rate  constant  defined  in  2q  22, 

was  the  major  parameter  with  which  previous  investigators 
characterized  the  combustion  process  of  solid  fuel  particles. 
This  parameter  was  calculated  in  Appendix  I to  be 
^diff  = 2 ft  A*diff  do2  at  T = 2320°K.  For  particles  with  an 
average  diameter  of  one  micron,  was  found  to  be  equal 

to  9^9  /*-2/sec  for  a diffusion  limited  reaction.  Particle 
burning  .times  measured  by  Macek  and  McLain  (Refs  9 and  11) 
suggested  an  evaporation  rate  constant  of  2,3  x 10'*  to 
9 x 1C1*  /A?/sec  or  twenty  to  one  hundred  times  greater  than 
that  given  in  Appendix  I.  Particle  agglomerations  in  the 
present  investigation  could  account  for  much  of  this 
discrepancy  and  the  average  combustion  temperature  could  also 
have  some  effect.  Assuming  that  the  differences  in  the 
combustion  temperature  and  the  agglomeration  effects  were 


DS/ME/74-2 

small,  the  agreement  In  the  evaporation  rate  constants  would 
be  much  better  if  the  reaction  rate  process  was  limited  by 
the  chemical  reaction  and  not  by  diffusion.  This  would 
yield  evaporation  rate  constants  from  400  to  2000  /*/sec 
from  Macek's  and  McLain's  data  and  4?4  )U/sec  from  the 
present  investigation.  There  is  some  evidence  that  this  may 
be  the  case.  Macek's  data  at  2280°K  and  2430°K  gas  tempera- 
tures for  the  35  and  44  micron  particles  correlates  best  for 
the  linear  burning  law.  Sheer  (Ref  26)  measured  the 
reaction  rate  x>f  99$  pure  boron  powder  at  425°C  to  495°C  and 
reported  ah  activation  energy  of  40  Kcal/mole  for  that 
reaction.  This  reaction  temperature  was  low  enough  that  the 
reaction  rate  should  have  been  limited  by  chemical  reaction 
rather  than  diffusion.  This  was  the  same  value  for  the 
activation  energy  deduced  in  the  present  investigation. 

Since  there  is  no  significant  difference  in  the  analysis 
between  the  chemical  limited  and  the  diffusion  limited  rate 
processes  for  a one  particle  size  distribution,  there  was 
no  way  to  resolve  this  question  in  this  investigation. 

Because  the  diffusion  process  was  the  concensus  of  the 
previous  investigators,  it  was  used  in  the  present  analysis. 
Diagnosis  of  Past  Combustion  Problems 

All  attempts  in  the  past  to  use  boron  as  a fuel  additive 
have  had  difficulties  in  making  the  boron  burn  within  the 
combustor.  These  problems  have  been  the  difficulty  in 
Igniting  the  boron  and  slow  boron  reaction  rates  (Ref  5). 

The  results  of  this  investigation  established  residence  times 
for  stable  combustion  and  complete  oxygen  consumption  to  be 
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approximately  two  milliseconds  and  four  milliseconds, 
respectively. 

The  improvements  observed  in  this  investigation  point 
out  the  two  main  difficulties  that  have  plagued  past 
attempts  to  obtain  high  combustion  efficiency  at 
atmospheric  pressure— particle  agglomerations  and  energy 
losses  from  the  flame.  The  agglomeration  problem  was 
solved  by  the  dry  cloud  of  boron  particles  from  the  particle 
mill.  The  dispersed  particles  decreased  the  particle 
agglomerations.  The  Insulated  combustors  with  the  boron 
and  its  oxide  deposits  on  the  combustor  walls  reduced  the 
heat  losses  from  the  flame  and  aided  in  the  rapid  heat  up 
to  ignition  of  the  suspended  fuel  particles.  The  reduction 
of  energy  loss  from  the  flame  enabled  the  flame  to  burn  at 
a higher  tempera tu"  This  caused  the  reaction  to  proceed 

faster  and  eliminated  the  inhibiting  oxide  coating  on  the 
fuel  .particles. 

The  method  of  preparing  the  fuel  aided  in  providing 
cleaner  fuel  particles.  The  boron  oxide  is  soluble  in  water 
and  was  washed  away  when  the  boron  powder  was  pressed  into 
solid  pellets.  The  closely  packed  pellets  did  not  oxidize 
as  rapidly  as  a dry  packed  powder  would  have.  The  wire 
brush  particle  mill  stripped  these  cleansed  particles  from 
the  pellet  surface  and  delivered  them  to  the  combustor 
before  they  had  time  to  oxidize. 
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V.  CONCLUSIONS 

This  research  was  conducted  to  aid  the  understanding 
of  multiple  particle  combustion  of  boron.  The  specific 
objective  was  to  deduce  an  analytical  combustion  model  for 
steady  state,  self-sustained  combustion  of  an  air/boron 
particle  mixture  and  to  demonstrate  its  validity  experi- 
mentally, A combustion  model  was  deduced  that  accounted 
fors  (a)  the  energy  mechanisms  required  to  heat  the  gas- 
particle  mixture  to  the  ignition  temperature,  (b)  the 
energy  loss  mechanism  from  the  flame  zone  by  unconsumed  fuel 
particles  and  the  combustion  products,  and  (c)  the  tempera- 
ture dependent  particle  reaction  rate. 

The  steady  state,  self-sustained  combustion  of  boron  in 
air  was  established  and  investigated  to  evaluate  the  param- 
eters identified  in  the  analysis.  The  experimental  combus- 
tion process,  operating  in  the  stoichiometric  to  fuel  rich 
regime,  achieved  high  values  of  oxygen  depletion  and  energy 
release. 

From  this  investigation,  it  was  concluded  that* 

1,  The  analytical  model  for  one-dimensional  flow 
agreed  with  the  measured  characteristics  of  temperature  and 

oxygen  depletion  for  particulate  boron  combustion  in  constant  j 

area  combustors  at  combustion  pressure  near  atmospheric,  ; 

2,  The  combustion  system  design  techniques  of  dry 

powdered  boron  and  insulated  combustors  used  in  this  study  | 

fl 

resulted  in  rapid  burning  and  high  combustion  efficiency,  ij 

¥ 

Complete  depletion  of  oxygen  within  four  to  five  milliseconds  \ 

was  measured  for  equivalence  ratios  ^ between  one  and  two,  $ 

si 

i 
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3.  The  primary  combustion  product  of  boron  and  air  was 
B2O3.  This  was  indicated  by  measured  combustion  temperatures 
near  the  adiabatic  flame  temperature  for  BgO^  and  above  the 
temperatures  possible  for  a chemical  equilibrium  process. 

4.  The  reaction  rate  for  the  particle  burning  process 
was  temperature  dependent  and  followed  an  Arrhenius  equation 
with  activation  energy  of  40  Kcal/mole  and  pre-exponential 
constant  of  9.04  x 10“3  moles  of  air/om-sec. 

5.  The  evaporation  rate  constant  for  the  particle 
correlated,  best  for  chemical  limited  rather  than  for 
diffusion  limited  combustion.  The  value  for  the  chemical 
limited  reaction  was  calculated  to  be  474  /t/sec.  For 
diffusion  limited  combustion  this  constant  was  calculated  to 
be  949  /42/ sec. 

6.  When  the  gas -particle  mixture  became  luminous  from 
combustion  heating,  the  absorption  coefficient  K was 
increased  ten  fold. 

7.  The  highest  flame  luminosity  occurred  at  flame 
temperatures  just  above  the  minimum  necessary  for  self- 
sustained  flames  (2100  to  2200°K). 

8. . The  calculated  absorption  coefficient  for  the  flame 
was  proportional  to.  the  absorption  coefficient  measured  at 
the  0.518  P-  wave  length. 

9.  The  flame  absorption  coefficient,  measured  at 
0.518/4,  was  proportional  to  the  theoretical  reaction  rate 
of  the  combustion  process  which  suggested  the  primary 
absorbing  species  may  be  a short-lived  intermediate  species. 
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10.  The  absorption-emission  technique  for  deducing 
high  temperatures  provided  an  accurate  method  for  the 
measurement  of  luminous  multiphase  flame  temperatures.  The 
oscilloscope  traces  of  the  flame  emissions  had  large  fluctua- 
tions in  intensity;  but,  the  maximum  level  of  the  signal 
represented  the  temperature  of  the  flame. 

11.  Self-sustained  boron/air  combustion  in  constant 
area  combustors  was  not  possible  without  some  of  the  boron 
adhering  to  and  burning  on  the  combustor  walls.  However,  at 
very  low  air. mass  fluxes  (on  the  order  of  0.1  gm/cm2-sec)  it 
was  possible  to  propagate  the  flame  into  the  cold  gas-particle 
mixture. 

12.  Thickening  of  the  wall  deposit  at  the  combustor 
entrance  occurred  for  air  mass  fluxes  less  than 
1*5  gm/cm2-sec;  this  resulted  in  decreased  combustor  volume. 
Air  mass  fluxes  above  1.5  gm/cm2-sec  produced  uniform  wall 
deposits. 

13.  Past  combustion  problems,  the  difficulty  to  ignite 

the  boron  particles  and  the  slow  combustion  reaction  rates, 
were  resolved  with  this  combustion  system.  These  problems 
were  primarily  due  to  particle  agglomerations  and  energy 
losses  from  the  flame.  The  agglomeration  problem  was 
solved  by  the  dry  cloud  of  dispersed  particles  from  the 
particle  mill  and  the  combustors  constructed  of  high  thermal 
resistant' materials  with  the  boron  and  boron  oxides  on  the 
combustor  walls  reduced  the  energy  losses  from  the  flame 
and  aided  in  the  rapid  heat  up  to  ignition  of  the  suspended 
particles.  \ 
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RECOMMENDATIONS  FOR  FURTHER  INVESTIGATIONS 

As  in  any  investigation,  there  were  almost  as  many 
questions  raised  as  there  were  answered.  The  following 
topics  are  among  those  questions  that  have  arisen  in  the 
course  of  this  study. 

1.  Determination  of  the  pressure  dependency  of  the 
boron/air  system.  The  present  investigation  was  conducted 
at  atmospheric  pressure  but  most  combustion  systems  of 
interest  are  at  either  subatmospheric  pressure  for  high 
altitude  flight  or  at  elevated  pressures  for  sea  level,  high 
Mach  flight  speeds.  The  boiling  point  of  boron  trioxide 
would  be  changed  and  could  change  the  reaction  limiting 
phenomena.  Increased  air  density  would  cause  a correspond- 
ing increase  in  the  fuel  loading  per  unit  volume,  this 
could  shift  the  optical  limit  toward  the  thick  region. 

Thus,  this  combustion  system  could  exhibit  very  different 
properties  during  high  pressure  combustion. 

2.  Evaluate  larger  combustor  geometries  to  determine 
the  changes  in  the  effects  of  the  wall  deposits,  wall 
reactions  and  flame  emissions.  An  increase  in  the  combustor 
diameter  affects  the  ratio  of  the  combustor  volume  to  surface 
area,  which  could  affect  the  wall  reaction  to  flame  emission 
ratio  for  the  energy  required  to  heat  up  the  cold  particle 
cloud.  The  Reynolds  number  change  for  the  same  air  mass 
flux  could  change  the  character  of  the  wall  deposits  and 
reactions.  The  two-dimensional  effects  caused  by  the  larger 
diameter  particle  cloud  could  completely  alter  the  combustor 
length  parameter. 
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3.  Further  investigation  of  the  absorption-emission 

technique  for  measuring  the  temperature  of  luminous  flames 
containing  liquids  and  solids  is  needed  to  evaluate  the 
accuracy  of  the  system  in  these  environments.  The  assump- 
tions made  in  the  present  study  were:  (a)  thermal  equilib- 

rium between  the  particles  (liquid  or  solid)  and  the  gases 
in  the  flow,  (b)  translational  and  molecular  temperature 
equilibrium  of  the  gases,  (c)  the  temperature  measurements 
taken  at  the  gaseous  emission  peak  were  not  affected  by  the 
particle  .continuum  emissions,  and  (d)  scattering  and 
reflection  extinction  and/or  enhancement  mechanisms  were 
neglected.  The  comparison  of  flame  emissions  and  calibrated 
source  attenuation  phenomena  at  different  wave  lengths  may 
aid  in  determining  the  validity  of  these  assumptions.  Of 
particular  interest  would  be  the  comparison  of  the  data  at 
wave  lengths  where  the  gaseous  species  are  highly  absorbing 
to  wave,  lengths  where  the  gases  are  transparent.  These 
measurements,  in  conjunction  with  alternate  measurement 
techniques  would  be  necessary. 

4,  Evaluation  of  minimum  length  dump  combustors  or 
other  recirculating  flow  geometries  on  reaction  efficiency. 
This  system  may  require  two-dimensional  and  three-dimensional 
analyses  of  multiphase  reacting  flows.  The  spherical  parti- 
cle assumption  may  be  applicable  to  thermal  radiation 
absorption  but  the  irregular  shapes  of  the  particles  should 
be  considered  for  the  solid  particle  species  dispersion  and 
the  momentum  transfer  between  the  gases  and  solids.  A 
simpler  analysis  could  be  to  correlate  the  combustor 
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geometries,  the  gas  flow  properties  and  the  solid  loading  <f> 
to  the  minimum  combustor  length  necessary  for  high  combus- 
tion efficiency. 

5.  Investigate  the  B02  concentration  in  boron  flames 
to  determine  its  relationship  to  the  reaction  rate  and 
ignition  delay  txme.  The  single  particle  investigations 
have  correlated  the  visible  emissions  from  the  particle  to 
the  minimum  ignition  temperature,  ignition  delay  and  to  the 
time  necessary  to  completely  consume  the  particle.  Kaskin 
and  John  (Ref  J24)  have  Identified  these  visible  emissions 
to  eminate'  from  the  B02  molecule  and  net  from  the  final . 
combustion  product  B20-j.  The  present  Investigation  has 
shown  the  absorption  coefficient  of  these  emissions  to  be 
proportional  to  the  boron/air  reaction  rate.  An  investi- 
gation into  these  phenomena  could  yield  understanding  of 
the  order  of  the  reaction  rate  kinetics,  the  temperature  and 
concentration  of  B02  required  to  give  visible  emissions  and 
the  phase  relationships  between  B02  and  B20^  in  nonchemical 
equilibrium  reacting  systems. 

6.  Determine  the  dominant  mode  for  limiting  the  combus- 
tion of .multiparticle  boron/gas  mixtures.  This  investiga- 
tion would  require  uniformly  distributed  particle  clouds 

of  different  sized  fuel  particles  or  particle  distributions 
with  different  average  diameters.  Previous  investigators 
have  assumed  a diffusion  limited  process  for  the  combustion 
of  single  particles  while  the  chemical  limited  process 
correlates  the  present  data  best.  Talley  (Ref  15) 
characterized  the  combustion  process  for  large  particles  of 
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boron.  He  found  that  below  the  boiling  point  of  BgO^  the 


rapid  reaction  of  the  boron  was  limited  by  the  diffusion  of 
the  oxide  vapor  away  from  the  reacting  surface.  Once  the 
boron  was  free  of  condensed  oxides  on  its  surface,  the 
combustion  limiting  process  would  be  either  chemical  or 
diffusion  of  the  oxidants  to  the  reacting  surface.  At  much 
higher  temperatures,  the  diffusion  of  boron  vapor  away  from 
the  surface  would  be  a limiting  process. 
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APPENDIX  A 


CHEMICAL  EQUATION  AND  DEFINITION  OF  CHEMICAL  BATIOS 
The  molar  chemical  equation  assuming  no  nitrogen  reac- 
tions and  air  composed  of  oxygen  and  nitrogen  only  is 


a,  & + *2  (Ox+Xo  Nt  ) “**  a<f  6t  0^  f B t OfO^  CL&  fy,  (A-l) 


The  stoichiometric  equation  would  then  be 


(«,*%) 8 +&-«#Vq  +^»  NP  ——  o.ij (<tz-a.JX,Nt  u-z) 


B.  i 


From  these  equations,  the  following  terms  may  be  defined: 


(840-  O-Mb 


(A-3) 


(B/A>,=  (*'- 


(a-i>) 


(A-5) 


>?  = [(  rr)0l~  >«/>]/(%)„  = ^,-Oj)/a,  (a-6) 
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/ - ^/(m.J'  = («2-  «4  )/a2  (A-?> 

Combining  Eqs  A-5,  A-6  and  A-7,  one  obtains 


V s / - n^/cm^  - <j>y\  (A_8) 

Since  BgO^  is  the  only  combustion  product  considered,  a^j 
is  equal  to  zero  for  all  i,  ;}  except  1=2  and  J = 3.  For 
this  case,,  (a^.  - a^)/(a2  - a^)  = 4/3*  The  fuel  used  in 
this  investigation  had  10#  impurities  by  weight;  thus,  the 
actual  fuel-to-air  stoichiometric  ratio  was  110#  of  the 
boron/air  ratio.  The  boron/air  ratio  by  weight  for  BgO^ 
as  the  only  product  is  (B/A)s  = 0.1046,  then  (F/A)s  = 0.115. 


1 


DS./ME/74-2 

APPENDIX  B 

GAS  PHASE  VOLUMETRIC  AND*  K/>SS  CHANGES 
The  change  in  the  mass  and  molecular  weight  of  the 
boron/air  reaction  oan  be  calculated  using  the  chemical 
equation  for  an  incomplete  reaction  quenched  at  some  time 
t after  the  start  of  the  reaction. 


M0  B + (N0^  Qg + (N„t)  S/2  ->  A/^  Bj,03  + N0  B + 0/(^\  < b-1  > 


Where  Nm  is  the  number  of  moles  of  species  m. 

For  conservation  of  the  elementary  chemical  species, 
the  right  side  of  Eq  B-l  is  subject  to  the  following 


constraints: 


= (Ne\ 

3/z%03+  *oz-Wo2l 


(B-2) 


Using  the  definitions  for  the  combustion  efficiency 
and  oxygen  depletion  parameter  ^ , introduced  in  Appendix  A, 
the  first  relation  of  Eq  B-2,  the  total  number  of  moles  on 
the  rigljt  side  of  Eq  B-l  can  be  expressed  in  terms  of  the 
initial  molar  conditions 

Hr  --  Ki  l ♦ Kl  m - W.  * (Sl 

= Kr  X * - * ("el  - < *>  + ' WflJL  <B-3) 
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Where 

M«l  = ^Ajl  + Wki 

and 

Va=  **fSi 

Due  to  its  high  boiling  point,  3930°K,  boron  can  be 
assumed  to  remain  in  the  solid  state  with  negligible  partial 
pressure  throughout  the  entire  reaction  process.  However, 
B2°3*  wili°k  bolls  at  2320°K,  would  have  a significant 
partial  pressure  and  could  be  completely  vaporized  at  high 
reaction  .efficiency.  Thus,  the  number  of  moles  in  the  gas 
phase  in  Eq  B-3  is 

1%  = (N a,rl~  O-ft^Wct). 

= [l-(/-%P%J(C*,).fHrl  (B-4) 

Where  P_  was  the  partial  pressure  of  BoOo  and  ranged 

B2°3  J 

from  zero  to  one  and  (C 0 )Q  was  the  mole  fraction  of  oxygen 

in  air  with  a value  of  0.21. 

At  low  values  of  the  oxygen  depletion  parameter  'Y'  , the 

partial  pressure  of  BgO^  would  be  small.  At  high  values  of 

reaction  efficiency  ( ^^1),  the  partial  pressure  would 

approach  unity.  Thus,  Eq  B-4  can  be  written 


<Jr.=  0- 6.llf)(Nair\  ;f<<' 

= (l-0°TP)(Na,r\  ;1 (-*/ 


(B-5) 


\ 
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The  change  of  mass  within  the  gas  phase  can  be  calculat- 
ed by  multiplying  the  number  of  moles  by  its  molecular  weight* 
From  Eq  B-4,  this  becomes 


= ( /"  O.ti'fX  >»«,„>„  ; V«  I (B-6) 

= (1+  o.i&WC  j y a) 


Where  (02/air)o  is  the  mass  fraction  of  oxygen  in  air  with  a 
value  of  0.23  and  is  the  molecular  weight  of  species  m. 

The  changes  in  the  equation-of -state  for  the  gaseous 
elements  can  be  determined  using  Eq  B-4 


(B-7) 


Frpm  this  analysis,  the  volumetric  and  mass  character- 
istics of  the  gaseous  elements  do  not  vary  significantly 
from  that  of  air.  For  gaseous  the  maximum  changes 

were  seen  to  be  7#  with  respect  to  volume,  10#  in  mass  with 
an  Increase  in  the  molecular  weight  of  12#.  For  no  gaseous  | 

B203,  the  mass  and  volumetric  changes  would  be  at  a maximum;  jj 

however,  there  would  be  almost  no  change  in  the  molecular  1 

weight.  \ 1 

t _ — 
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APPENDIX  C 

DEVELOPMENT  OF  THE  CONSERVATION  EQUATIONS 


The  conservation  equations  for  a flowing  system  were 
developed  using  the  Reynolds  Transport  Theorem  (Ref  31 )• 

The  Reynolds  Transport  Theorem  states  that  for  a quantity  F 
that  can  be  described  per  unit  mass  of  tne  substance  within 
a closed  volume,  the  change  of  F with  respect  to  time  within 
the  volume  is  equal  to  the  change  of  F at  each  point  within 
the  volume  plus  the  flux  of  F across  the  boundary  surfaces  S 
containing  the.  volume . Stated  in  vector  form,  the  Reynolds 
Transport  Theorem  is 


For  steady  flow  conditions,  there  is  no  change  in  any 
quantity  within  a volume  fixed  in  space,  thus,  for  steady 
flow,  Eq  C-l  is 

A JfJ  S’?  ^ <°-2> 

V 5 

Invoking  the  Divergence  Theorem,  the  surface  integral  may  be 
stated  as  a volume  Integral 


(C-3) 


In  a multispecies,  multiphase  system,  each  phase  or 
species  within  the  phase  may  have  a different  velocity. 

Thus,  for  a flow  system  of  this  type,  the  flow  properties  in 
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Eq  C-3  must  be  expressed  in  terms  of  each  species  m and 
phase  j 

(o'4) 


The  Conservation  of  Mass 

Since  there  are  no  sources  or  sinks  of  mass  even  In  a 
reacting  system,  the  mass  Is  conserved  and  the  Integrals  of 
Eq  C-4  are  identically  zero  for  P = 1.  Then 


J/zz^.rs  = o 


(c-5) 


For  one -dimensional  flow,  the  only  velocity  component 
is  the  axial  velocity  u and  the  flow  properties  are  constant 
across  any  cross  section  perpendicular  to  the  flow,  then 
Eq  C-5  can  be  integrated  to  yield 

or 

and  for.  a constant  area  duct 


= Const  t 


(C-6) 


When'  the  velocities  of  all  species  and  phases  are 
assumed  to  be  identical,  there  is  no  slip  between  the  solid 
or  liquid  particles  and  the  carrier  gas.  This  implies  that 
there  was  conservation  of  the  elementary  molecular  or  atomic 
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species  as  well  as  conservation  of  mass.  Thus,  there  was  no 
change  In  the  fuel/alr  ratio  within  any  given  volume  fixed 
In  the  combustor  space.  For  a reacting  flow,  the  Individual 
chemical  species  are  not  conserved  but  the  sum  of  all 
species  m Is  conserved.  Thus,  Eq  C-6  can  be  written 


- Consent 


(C-7) 


Initially,  before  any  reactions  take  place,  the  mass 
of  the  gas.-par.ticle  mixture  is  solid  boron  particles  and 
air 


(0-8) 


= 0 + FA)(/Sr\H. 


Since  the  fuel/alr  ratio  remains  constant  within  the  combus- 
tor space,  Eq  C-7  can  be  expressed  In  terms  of  the  air 
before  the  reaction  had  started.  •* 

(C ).K.=  v I r /=}/  (I* r/A)  = /f„u  = C0«1«1  (0-9) 

The  Momentum  Equation 

The  rate  of  change  of  momentum  of  a flow  system  is 
equal  to  the  summation  of  the  forces  acting  upon  It.  For 
a horizontal  flow  system,  the  gravitational  forces  can  be 
neglected  and  neglecting  the  viscous  forces,  only  the 
pressure  forces  P acting  on  the  surface  of  the  volume  must 

be  considered.  The  momentum  per  unit  mass  is  the  velocity, 

#• 
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thus,  Eq  C-4  when  F = v is 


■ If 2 If. 


w>  0 


(C-10) 


Where  only  the  gas  phase  contributes  to  the  pressure  tensor. 

Using  the  Divergence  Theorem  and  rearranging,  Eq  C-10 
becomes 


v - * ’ r 

Assuming  that  all  properties  within  the  Integral  are 
valid  at  each  point  within  the  integral  and  noting  from 


representation  of  the  momentum  equation  C-ll  can  be  written 


° (0'1: 

Where  V'P  = -VP 

For  one-dimensional  flow  and  no  slip  flow  between  the 
phase  states,  Eq  C-12  becomes 


(C-13) 


Equation  C-13  must  be  summed  over  all  phases  not  just 
the  gaseous  one.  To  take  the  solid  and  liquid  phases 
separately  would  imply  that  these  velocities  are  zero  since 
they  do  not  contribute  to  the  total  pressure.  To  maintain 
the  velocity  of  the  particles  approximately  equal  to  that  of 


112 


^77  ^^e^^SWS^npTSJP^SP1®^*® 


nmMMRMR 


DS/ME/74-2 

the  gases  will  require  energy  extraction  from  the  gases* 

Thus,  the  solids  will  contribute  to  the  pressure  gradient 
in  multiphase  flows. 

The  Energy  Equation 

The  rate  of  charge  of  energy  within  a volume  is  equal 
to  the  energy  source  or  sink  caused  by  chemical  reaction 
^chem*  worlc  due  fco  pressure,  plus  the  heat  flow  into  the 
volume  and  minus  the  heat  flow  out  of  the  volume.  For  each 
phase  j of  a flow  system,  the  energy  per  unit  mass  of  species 
m is  (em^  f Y^/2)  and  since  the  heat  flux  Inward  is  positive, 
the  integral  energy  equation  can  be  written 


U{fl\w4>w- 


JV 

(C-14) 


Where  is  the  internal  energy  per  unit  mass  of  species  m 
in  phase  J and  q was  the  combined  heat  fluxes  shown  in  Fig  3» 
Using  the  Divergence  Theorem  to  express  the  surface 
integrals  as  volume  integrals,  Eq  C-14  can  be  written 


(C-l 5) 


Assumption  (b)  of  the  flow  model  states  that  the  number 
of  solid  particles  was  constant  in  any  given  volume  of  the 
gas-particle  mixture,  thus,  Eq  C-15  is  valid  for  every  small 
volume  as  well  as  the  integral  volume.  For  one-dimensional 
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flow  and  the  velocity  and  temperature  independent  of  the 
system  phase  states,  the  differential  element  of  Eq  C-15 
becomes 


££ Vi'v#*) + = ° (0-16> 

Expanding  the  pressure  work  term  in  Eq  C-16 

r ^ + *zz  i&i-  (0-17) 

The  derivative  of  the  velocity  can  be  obtained  by 
differentiating  Eq  C-9.  Substituting  Eqs  C-9  and  C-13  into 
C-17,  one  obtains 

z J(r>u) = -(i  e,)  * /&  - r (o-is) 

^ </*'  <7  ^ /Sr  ^ * f * M 

Equation  C-l?  can  be  expanded  further  using  the 
equation-of -state  for  air  and  rearranging  the  last  term. 


Substituting  Eq  B-7  into  Eq  C-19,  the  pressure  work 
term  may  be  expressed  as 

l &(&#>--  - f/-3 V §T 


-TT & uJ*% 

* jx 


(C-20) 


-v.  a 
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Substituting  Eq  C-20  into  Eq  G-l6,  collecting  terms  and 

using  the  definitions  (cv)s  = (op)s  = cs,  (cp)g  = (cv)g  + Rg 

and  demj/dx  = cy  dT/dx,  Eq  C-l6  can  be  simplified  to 
J m;| 
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APPENDIX  D 

DETERMINATION  OP  THE  VOLUMETRIC  SPECIFIC  HEAT 
The  total  volumetric  specific  heat,  is  a 

IW  ^ 9 O' 

measure  of  the  energy  stored  within  the  gas-particle  mixture 
per  degree  of  temperature  per  unit  volume.  Even  though, 
from  the  analysis  in  Appendix  B,  the  gas  phase  volume  and 
molecular  weight  did  not  vary  significantly  as  the  boron/ 
oxygen  reaction  proceeded,  the  specific  heat  capacity  of 
the  gas-particle  mixture  must  be  evaluated  to  determine 
whether  the  change  due  to  the  reaction  product,  BgO^,  was 
significant.  Since  the  mass  within  the  control  volume  for 
steady  state  condition  was  conserved,  the  atomic  species  was 
also  conserved  and  the  phase  states  can  be  expanded  into  the 
chemical  species  m.  Expanding  the  total  volumetric  specific 
heat  in  terms  of  the  mass  of  the  chemical  species  per  unit 
volume , 


The  fuel  mass  contained  impurities  as  well  as  boron. 

The  impurities  consisted  of  various  metal  oxides.  Table  III 
of  Appendix  J;  however,  the  percentage  of  each  Impurity  was 
small.  So  to  simplify  the  analysis,  it  was  assumed  that  the 
heat  capacity  of  the  fuel  was  identical  to  pure  boron. 

Now 


(D-2) 
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and 

"w  Xf  (n">)  V (B'3) 


From  conservation  of  the  elemental  species,  the  mass  of 
BgO^  can  be  expressed  in  terms  of  the  change  in  the  mass  of 
boron  and  oxygen.  From  Eq.  B-2,  the  mass  of  BgO^  is 


(D-4) 


Using"  the  definition  for  the  oxygen  depletion  parameter 
^£/,  derived  in  Appendix  A,  Eq  D-Jf  can  be  written 


”Vr  2/5 * (N^ w«a 


= Vsf 


n 


(D-5) 


= lA5V{m0jc 


Since  Mg  o =69.6  and  M q Is  32  gm/mole. 

2 3 2 

Substituting  these  relations,  obtained  in  Eqs  D-2,  D-3 
and  D-5,  into  Eq  D-l 

- w0,r 


Dividing  by  the  volumetric  specific  heat  of  air,  the  ratio  of 
specific  heats  © is  defined  as 
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®i^P,/%AW  (D- 

= /+r%)«a  - «/o  >1  p +(yA)  -%.  1 

€ir  £/»*  *“  ^ £//•“* 

Where  (02/A)  * {m0  )0/^mair^o* 

C* 

Using  the  definitions  (B/A)  * (B/A)s^>  and  y>  . 

derived  in  Appendix  A,  the  ratio  of  the  specific  heats  ©, 
Eq  D-6  may  be  written 


<2> 


+ v((W)[-4S 


(D-7) 


In  evaluating  Eq  D-7  numerically,  the  following  values  were 
used  (Ref  22): 


(P/A)s  * 0.115; 
(02/A)  = 0.23  ; 
Then 


CPair~  0,257  cal/®m“K; 

Cr,  =0.69  cal/gm-K; 
PB 


cp  = 0.353  cal/gm-K! 

cD  = 0.173  cal/gm-K. 

P°2 


© = 


/ +0.31  £ ■+'&  ( o.2i  [a si  - o.n. ?J  -Co./o«X<>.ei}^o.zs7 
l + 0.il<l>  +0.02IV 


The  third  term  was  the  change  of  heat  capacity  per  unit 
volume  of  the  mixture  due  to  combustion.  Since  it  was  small 
compared  to  the  rest  of  the  relation,  it  was  neglected. 

Thus,  the  value  for  ©,  used  In  the  combustion  model  was 


©=  l+0.5l4> 


(D-8) 
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APPENDIX  E 

RADIANT  ENERGY  EQUATIONS  OP  TRANSFER  (REF  20) 

The  gas-particle  mixture  was  assumed  to  be  a flow  of 
solid  particles  suspended  in  the  gas  and  uniformly  distrib- 
uted throughout  the  gas  volume.  Before  vigorous  combustion 
of  the  particles  took  place,  the  gas  was  transparent  to 
radiation  and  only  the  solid  particles  absorbed,  scattered 
and  emitted  radiation.  During  vigorous  combustion,  the  solid 
fuel  particles,  the  solid  combustion  products,  as  well  as  the 
gaseous  combustion  products  could  participate  in  the  radiant 
transfers.-  Thus,  a cloud  of  solid  particles  suspended  in  a 
gas  volume  would  behave  like  an  absorbing,  emitting,  scatter- 
ing gaseous  media. 

During  the  heat  up  phase,  the  gas-particle  mixture 
would  be  bounded  on  the  upstream  end  by  a black  cavity 
radiating  at  the  ambient  temperature  T0  and  on  the  downstream 
end  by  a radiating  combustion  zone  that  was  effectively  a 
black  cavity  radiating  at  an  average  flame  temperature  Tp, 

Pig  1.  Once  the  combustion  zone  was  reached,  the  major  heat 
flux  from  an  elemental  layer  of  the  mixture  was  upstream 
toward  the  heat  up  zone  and  the  radiant  transfer  downstream 
into  the  exhaust  plume  was  neglected. 

To  evaluate  the  radiant  heat  fluxes  in  an  absorbing, 
emitting  and  scattering  media,  consider  an  infinitesimal 
layer  of  the  media  lying  between  the  distances  x and  x + dx 
from  a cavity  radiating  at  T0.  This  media  will  be  assumed 
to  be  a homogeneous,  gray  media  with  isotropic  coefficients 
of  scattering  s,  emission/'  and  absorption  k. 
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The  radiant  energy  transfer  balance  can  be  obtained 
from  a beam  of  energy  I of  solid  angle  du  as  it  passes 
through  the  layt  of  the  media  at  some  angle  © from  the 
normal  direction.  Pig  27.  The  energy  within  the  solid  angle 
was  changed  by  a differential  dl  as  it  passed  through  the 
distance  dx/Coso  • The  change  was  due  to  the  amount  of 
energy  in  I absorbed  in  the  layer  of  the  media,  kl;  the 
amount  of  energy  in  I scattered  out  of  the  solid  angle  dw  , 
si;  the  energy  from  all  other  energy  beams  scattered  into 
the  solid,  angle,  s^f  J ; and  the  energy  emitted  from  the 
layer  of  media  along  the  solid  angle,  This  energy 

change  can  be  written 

Cose  dl  = -Jtl  -si  + v<r  T4  + s Joj 
Ax  tt  Wii  Jia 


4V 


= -Jl  7 4 -h  S f JJu> 

~ 4TT  Jfjt 


(E-l) 


TT 


Where  J.=  k + s and  is  called  the  extinction  coefficient. 
Equation  E-l  can  be  rearranged  to 


Cose  J2  + I = / IJu)  (E-2) 

j J7T  347 Un 

» 

Changing  the  independent  variable  to  the  optical  length  ^ 
where 

J Jx 

and 

r=  <E“3) 


\ 
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Multiplying  Eq  E-2  by  the  integrating  factor  6 


'*/CoS9 


one  can  obtain 


z£ir 

IV 


4nT4fi 


#*] 


'Z/COSQ 

€ 

COS& 


(B-4) 


The  energy  transferred  along  beam  © can  go  both  ways.  To 
differentiate  which  way  along  the  beam  the  intensity  was 
going,  the.,  convention  I+  denotes  energy  moving  in  a positive 
x direction  and  1“  for  the  negative  direction.  Thus,  for 
the  solid  angle  = J % , I+  integrated  over 

was  the  total  positive  direction  flux  and  1“  integrated  from 
I < 0 < 7T  was  the  flux  in  the  negative  direction. 

Integrating  I+  from  x = 0 to  x yields 


e / [f  T V fig  (B-5a> 


CosO 


and  1“  from  xp  to  x 


(Tf-V/cos*  i fr^<r^4,  s frj 

I =Ire  -wj  [7  7 +4T,P“ J H <«-»> 


The  total  longitudinal  flux  of  radiant  energy  can  be 
obtained  from  Eqs  E-5a  and  E-5b  by  summing  the  components 
of  the  radiant  intensities  over  the  appropriate  solid  angle 
hemispheres.  Since  the  media  was  assumed  to  be  homogeneous, 
there  was  no  dependence  on  § , the  radiant  energy  flux  qr 
in  the  x direction  can  be  written 
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r~-M 


2*  % i*  IT 

- f Jl  CoSO  Smod&dG * J f 1 Cos BStnodQ  </5 

o % 


JX/e 


*/Z 


-!"Jl*ecoseM**»*l/[^T4+£ s /»e rf® ty  (E‘6) 

*6  o O 


* 


i&  r **  */* 


Differentiating  Eq  E-6,  one  obtains 


V ^ 


7T 


ji  • /*  ~^/coi  & f J^rV/cos9  ^ 

Snedo+znJ  IF^  s/nodo  -4j*t 


r m 

~ \ll  Ju + z/b'T%AlIJu]l 


Air 


M 

ll 

cose 


(E-7) 


7T 


r *n  ryt  c 


cose 


Smodoclyi 


The  value  for  the  energy  scattered  into  the  solid  angle 
can  also  be  obtained  from  Eqs  E-5a  and  E-5b 


mi 


IT/2 


fl  Jo-Z  n-/l„e  iT/eia*  *f  [t  T * J—  ^smtdtofy 


ATT 


(E-8) 


(rr-v)/co*o  fr u<r 
e S/aoete~2J  j 
t L 
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For  an  optically  thin  media,  that  is  for  <<  1, 

and  noting  that  IQ  and  Ip  erainate  from  nonreflecting  cavities 
I0  = cT,,4/#  and  Ip  = , Eqs  E-6,  E-7  and  E-8  can  be 

written 

%(t)  = $r&  ~ «r  7*-  r (E-6a) 

- ^ = z<rX *+  *r£4-  ^ - f/iX  (E-7a) 

and 

flJcjs  2<rT*+ZO‘l^4  (E-8a) 

*rr 

Substituting  Eq  E-8a  into  E-7a  and  using  the  definitions 
from  Eqs  E-l  and  E-3,  Eq  E-7a  can  be  written  as  the 
longitudinal  component  of  the  radiant  energy  divergence 


- = a i UT0\  z&<rT;*-  4vr  T*  (E-9) 

The  resulting  formulations  for  the  radiant  flux  and 
divergence  are  seen  no  longer  to  depend  on  the  scattering 
coefficient  or  the  integral  equation.  This  simplified  the 
energy  equation  significantly.  Also,  the  emittance  and 
absorption  coefficients  which  are  dependent  on  the  tempera- 
ture of  radiation  for  a nongray  media  are  dependent  only  on 
the  temperature  of  the  differential  layer  T(x)  for. a gray 
media,  k(T0)  ••=  k(Tp)  s k(T)  and  by  Kirchhoff's  law 
k(T)  = ^ (T).  The  absorption  coefficient  was  assumed  to  be 
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dependent  on  the  temperature  of  the  media  and  not  a constant 
by  the  formulation  of  the  definition  of  the  optical  distance 
r in  Eq.  E-3 

r - f A Jt  4 4x 
•» 
o 

\ 

\ 

5 

i 
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APPENDIX  P 

VELOCITY  AND  TEMPERATURE  RELATIONSHIPS 
FOR  THE  GAS-PARTICLE  MIXTURE 
To  determine  the  degree  of  error  Introduced  in  neglect- 
ing the  Mach  number  effect  on  the  velocity- temperature 
relationship  in  the  gas-particle  mixture,  the  relation  for 
the  velocity  as  a function  of  Mach  number  and  temperature 
Mas  derived  assuming  identical  velocity  and  temperature  for 
all  species  and  phase  states.  The  conservation  equations, 
the  equatiqn-tof -state  and  the  definition  for  the  Mach  number 
are  used. 

The  Continuity  Equation 

The  continuity  equation  for  the  gas-particle  mixture  is 
given  in  Eq  C-9  as 


<c-9> 

♦ 

Differentiating  Eq  C-9,  one  obtains 

J.  __  _ J_  J&r  (F-l) 

k dz  ~ ~ (%tr  Jx 

The  Momentum  Equation 

Substituting  Eq  C-9  into  Eq  C-13  and  expanding,  the 
momentum  of  the  mixture  can  be  written  as 


= 0+m£h,v.jh  = -IF  ifk* 


(F-2) 
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The  Gas  Equatlon-of -State  for  Air 

The  pressure  term  in  Eq  P-2  includes  only  the  gas  phase 
and  from  Eq  B-7  it  was  concluded  that  the  pressure  was 
approximated  by  that  of  air 


l Pmf*  tfr*T=  e, 


(P-3) 


i J&r  _ i J/mi»  , j_  Jr 

%ir  J * Jz  T dX 


(P-4) 


The  definition  for  the  Mach  number  M» 


m«=  ^Yrgr 


(P-5) 


i JCrMp  - _z  Jl  i Jr 

d%  V.  d*  ~ T dx 


(F-6) 


Substituting  the  continuity  and  momentum  equations 
P-1  and  F-2  into  the  differential  equation-of-stete,  Eq  F-4, 
and  utilizing  the  equation-of-state,  P-3,  a relation  involv- 
ing only  velocity  and  temperature  is  obtained 


[(-  0*FA)V/er]i^= 


(F-7) 


The  variables  in  Eq  P-7  cannot  be  separated;  however, 
the  bracket  on  the  left  hand  side  may  be  expressed  in  terms 
of  the  Mach  number.  Thus,  substituting  Eq  P-7  into  Eq  F-6, 
relations  involving  Mach  number  and  velocity  and  Mach  number 
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and  temperature  are  obtained 


i chi  _ MHiVJ#  . 

* dx  fr<i+9&wC]yt 


JL  _ 

yyil  na*rA)ti£ J fiftf 


(F-8) 


2-iTr  dMVjfr  _ f I _ 

r dz  [ mi+wmmI  J ytfi  ~ \fnt 


za+f/A> 

Hd4f/Aml 


(p_< 

i]  dx 


Expressed  in  logrithm  form,  Eqs  F-8  and  F-9  are 


ifl 


[/i  a-Ky^/nt] 


51“  . 1’ 


Thus, 


M 


yt f? 


[l  +^4  T 


y Mv 


— Coti5^«nt 


(F-10) 


— consent 

=vsn  ym3  1*t 


(F-ll) 


Dividing  Eq  F-ll  by  Eq  F-10,  one  can  obtain  the  temperature- 


velocity  relationship. 


[i+U+f/a)YH$]  T 


= Cons 


ianl 


± _ /+  fAfflft  _T 

A - /+a+m(mil  J 1 


(F-1.J) 


! 

* 

( 
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Evaluating  Eq  F-ll  at  the  initial  conditions,  one  obtains 


y Mm* 

&Ni\ 


l±  (i+ViWMh  % 

.i+a+wxYMi),  J ^ 


(P-lla) 


Equation  F-lla  can  be  rearranged  as  a quadratic  equation  in 
(1  + F/A)  ^ %2 


[(i+VaWmZ]' -i-  (z-  tu  (j+ftMMZV lx  1 *i-°  (p.13) 

(.  0*WffH}l  T ) 


The  variation  of  (1  + F/A)  ^ MN2  with  temperature  can 
now  be  evaluated  using  the  initial  conditions 
T0  = 300°K,  ^air  = F/A  = 0.115^  for  (f>-  1 and  2 and 

Tmax  = 3000°K.  Solving  Eq  F-13 


T 

300 

1000 

2000 

3000 

cuF/A)Yni 

0.017 

0.05 

0.115 

0.16 

l*(l*F/A\YMi 
t+  Of-F/MMnl  _ 

1 

1.03 

1.095 

1.14 

Thus,  it  was  seen  that  little  error  was  introduced  for 
exit  flame  temperatures  below  20 00°K;  however,  for  higher 
entry  Mach  numbers  and  higher  temperatures,  the  error  would 
be  significant  but  still  within  the  uncertainty  of  the  flow 
measurements • 
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APPENDIX  G 

TEMPFRATURE  MEASUREMENT 

The  measurement  o emperature  using  the  absorption- 
emission  technique  is  a noninterference  method  for  luminous 
flames  (Refs  27  and  28).  The  technique  requires  a monochro- 
matic measurement  of  the  emissive  power  of  the  flame  Ip^  ; 
the  attenuation  of  a calibrated  power  source  viewed  through 
the  flame  (absorption)  1^  ; and  a measurement  of  the 
emissive  power  of  the  calibrated  source  Irt_„  • Prom  these 
measurements t the  absorption  coefficient  of  the  flame  gases 

^F>  ls 


*FA= 


(G-l) 


Assuming  that  the  flame  system  is  in  radiative  equilibrium, 
the  spectral  emissivity  would  be  equal  to  the  spectral 
absorptivity.  Then  the  measured  emissive  powers  Ip^  and 


Ics^  may  be  described  by  Klrchhoff's  and  Planck’s  radiation 


laws, 


T _ C , 


(G-2) 


•csa  - 


*S(frp[V*TQi]-)) 


(G-3) 


Where  P^  is  the  spectral  constant  for  the  monochrometer. 

For  wave  lengths  less  than  If*.  and  temperatures  less 
T 

' >>  1.  And  by  using  the  ratio  of 


n V*  T 

than  3000°K,  e c 
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A tungsten  ribbon  filament  lamp.  General  Electric 
Company  type  30A/6V/T24  calibrated  by  the  National  Bureau 
of  Standards  for  brightness  temperature  at  0*665 ft.  as  a 
function  of  lamp  current  was  used  as  the  calibrated  source. 
The  power  supply  for  the  lamp  consisted  of  two  heavy  duty 
batteries  connected  in  series  to  nroduce  a twelve  volt 
output.  The  current  regulation  was  accomplished  by  dropping 
the  potential  across  a one  ohm  reostat  to  obtain  a constant 
current  through  the  lamp.  The  lamp  current  was  monitored  by 
measuring  the  current  across  a 50  mv,  50  ohm  shunt.  The 
current  was  displayed  on  a digital  electronic  potentiometer 
calibrated  with  the  shunt.  The  current  was  measured  in 
0.02  ampere  digits  to  within  0.025#. 

The  brightness  temperature  S*  was  converted  to  true 
temperature  by 


1 

1 

1 

'-4 


'cs  - 


5a. 


I+  fc/vO] 


Where  the  emissivlty  of  tungsten  for  the  wave  length  range 
from  0.45/<  to  0.68/tc  was  described  by  Larrabee  (Ref  29)  as 
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£ ft,T)  =0.465*-*  ISSSA  +0.26 IS X to  4 7-0.1$ 0S//6*at  ( G-6 ) 
Calibration  of  Spectrum 

The  Warner-Swasey  501  Scanning  Spectrometer,  with 
entrance  and  exit  slits  set  at  0.1  mm,  was  used  to  monitor 
the  spectral  range  from  0.4  to  0.65/*  . The  detector  was  a 
type  4473  photomultiplier  set  at  994  volts.  The  spectrum, 
as  displayed  on  the  oscilloscope,  was  calibrated  using  a 
mercury  arp  lamp  to  obtain  the  mercury  lines.  The  spectrom- 
eter was  then  monochroma tically  set  to  the  maximum  signal 
of  the  0.518/*  emission  peak  of  the  flame  spectrum. 
Measurements 

The  flame  temperature  was  measured  at  the  0.518/*  wave 
length.  The  lamp  emission  signal  was  chopped  at  one  second 
intervals  to  obtain  the  Ip^  and  measurements.  The 
lamp  measurement  Icg^was  taken  ;}ust  prior  to  the  test  run 
with  the  tungsten  lamp  current  and  oscilloscope  signal 
voltage  maintained  constant  throughout  the  test.  Schematic 
representations  of  the  signals  obtained  are  depicted  in 
Pigs  28. and  29. 
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APPENDIX  H 

FORTRAN  PROGRAM  THEORETICAL  COMPUTATIONS 


FOR  SOLUTION  OF  EQUATIONS  33  THROUGH  38 
(RUN  ON  CDC  6600) 


PROGRAM  ORMAND (INPUT , OUTPUT ) 

DIMENSION  Y(2),YD0T(2) 

COMMON  TO,TF,PHI 

READ  1000 , TO , TIMEF, DTIME 

PHI=0 * 

DO  100  1=1,2 
PRINT  2000 
PHI=PHI+1. 

TF=2900. 

*F(PHI.GT.l. ) TF=2400. 

TIME=0 . . 

.Y(l)=300. 

Y(2)=0. 

10  CONTINUE 

CALL  F(TIME,Y,YDOT) 

PRINT  2010 , TIME, Y( 1 ) ,YDOT( 1} ,Y( 2) , YD0T(2) 

CALL  RKDF ( TIME, Y. 2, DTIME) 

IF ( TIKE. LE. TIMEF)  GO  TO  10 
100  CONTINUE 
STOP 

1000  FORMAT(4E10.0) 

2000  FORMAT ( 1H1 , 15X, 4HTIME, 19X, 1HT, 12X, 8HDT/DTIME, 17X, 
13HPSI, 10X, 10HDPS I/DTIME ) 

2010  F0RMAT(5E20.5) 

END 

SUBROUTINE  F(TIME,Y, YDOT) 

DIMENSION  Y(2) , YD0T(2) 

COMMON  TO,TF,PKI 
T=Y( 1) 

PSI=Y(2) 

THETA=1.+.31#PHI 
F1=2.33*PHI*(TF/1000. )**4* 

1 ( 1 .+ ( TO/TF) **4-2* ( T/TF ) **4 ) /( THETA*T/1000 . ) 
F4=28.5*(1.+.58*(PHI-1.))*(3.-T/1000.)+2.33*PHI* 
1(1. -PSI/PHI ) ** . 666 

IF(T.GE.2100. ) Fl=-F4*( T/1000 • )**3/THETA 

F2=.014*PHI 

IF(T.GE.2000. ) F2=0. 

F3= ( 5650 . - 1850 . *EXP ( 19 . -44000 ./T) ) /THETA 
IF(T.GE.2320.)  F3=3800. /THETA 

IF(PSI.GE. .8)  F1=0. 

Gl=427 . *PHI* ( 1 . -PS I ) * ( 1 . -PSI/PHI ) ** . 33333 
G2=EXP(8. 62-20000. /T)/T 
YDOT( 1 ) =F1+F3*( F2+G1*G2 ) 

YDOT ( 2 ) =F2+G1*G2 
RETURN 

END  x 
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APPENDIX  I 
EVALUATION  OF  THE  CONSTANTS  A AND  fi 
The  reaction  rate  constant  and  the  Arrhenius  pre- 
exponential constant  A are  characteristics  of  the  reaction 
process  and  must  be  evaluated  from  experimental  data.  The 
reaction  rate  constant  /3  is  defined  for  a spherical 
particle  reaction  process  at  constant  temperature  and  un- 
limited oxygen  so  that  the  concentration  of  oxygen  CQ 
remains  constant  throughout  the  combustion  process.  The 
empirical  law  for  this  reaction  process  states  that  the 
rate  of  consumption  of  a fuel  particle  is  proportional  to 


the  exposed  surface  of  the  particle  and  the  concentration  of 
oxygen  surrounding  the  particle.  This  law,  for  a first  order 
reaction  rate  (Hef  30),  is  described  as 

_ JJUL  - t>tt  rz  Cox  Mr  JlSH.  (1-1) 

4h  4 ra 

Where  n is  the  number  of  particles  per  unit  volume  and  K(T) 
is  a constant  of  proportionality  ..^at  is  dependent  upon 
temperature.  The  value  of  a is  equal  to  zero  for  chemical 
limited,  reaction  and  equal  to  one  for  diffusion  limited 
reaction.  The  temperature  dependence  of  K(T)  has  been 
defined  in  exponential  form  by  Eq  25  as 

K(t)  = a e'E‘/RT  (25) 

Where  Ea  and  A are  defined  as  the  Arrhenius  activation  energy 
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and  pre -exponential  constant,  respectively,  and  Eq  25  as  the 
Arrhenius  equation. 

The  mass  per  unit  volume  of  the  fuel  particles  was 
expressed  in  terms  of  the  particle  volume,  its  liquid  or 
solid  density  and  the  number  of  particles  per  unit 
volume  in  Eq  7 


(7) 


Substituting  .Eq  7 into  Eq  1-1  and  expressing  the  relation 
as  a function  of  diameter,  one  obtains 

- A ( » 2L i.  ) 2*  n TT  4 i 

Integrating  Eq  1-2  yields 

i0  - i - 

Where  ft  - (<"*,)  2"'  Co*  Mf  K(T)/^ 

A similar  analysis  of  Eq  1-1  can  be  obtained  for  a 

finite  source  of  oxygen  and  a temperature  dependent  reaction 

process.  Dividing  Eq  1-2  by  the  initial  mass  per  unit  volume 

of  fuel  (mf)Q,  the  initial  concentration  of  oxygen  (CQ  )Q 

2 

and  K(T)  evaluated  at  some  reference  temperature  T^,  one 
obtains 


K(T)  (1-2) 


(1-3) 


A 

= / /sec 


\ 
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Where 


. A \(J/J]  - A'  (A  )2'a  Jk  K(t) 

V ft.,)  Kfr*> 

l'=  3[(*")2,‘*'(Ci!\Mf  KtWS^/fa  + OJi 

- 3 & 

(a+0  da*‘ 


(1-4) 


d-5) 


The  Arrhenius  pre-exponential  constant  A may  be  obtained 
from  Eq  1-5  by  expressing  K(Tk)  in  the  exponential  form 


/ - A'i(£"e£^ 

in, 


(1-6) 


The  experimental  value  of  A*  was  determined  from  Eq  37 

by  a graphical  analysis  of  the  data  in  conjunction  with  a 

correlation  of  the  data  using  the  solution  of  the  Fortran 

Program  listed  in  Appendix  H.  The  value  of  A'  that  agreed 

best  with  the  experimental  results  was  1423  sec~*  for  T 

1c 

equal  to  2320°K.  Evaluating  the  rate  constants  using  this 
value  and  an  initial  particle  diameter  of  one  micron, 
one  obtains 


/SJrt  = % Pi  £ - U 9 //sec 


(1-7) 


(1-8) 
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The  pre-exponential  constant  A can  be  evaluated  using  the 
following  values : 

Exp  [Ea/R232o]  = 5.55  X 103S  £B  = 2.34  gmB/cm3; 

(Cq  ) © = 0.21  N0  /Nair;  MB  = 10.81  gm/mole$ 

£ 4 

Ea/H  = 20,000»Ki  d0  = 10"1'  om. 


/4 


Wcur/C*»  -Sec. 


Noting  from  Appendix  A that  for  the  boron  reaction  in  which 

BgO^  is  the  only  product,  the  ratio  NB/NQ  = 4/3,  then 

2 

am f=  i-  04  K/o  3 Na,r/e  »•  - Sec  ( 1-9 ) 


and 

= lgl  y°'/cm3-  Stc.  (1-10) 


\ 
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APPENDIX  J 

TABLES  OP  PHYSICAL  PROPERTIES  FOR  METAL  FUELS 
AND  EXPERIMENTAL  RESULTS 
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afl 


A1 

933 

2767 

A12°3 

2315 

3250  - 

B 

2450 

3931 

b2°3 

723 

2320 

Be 

1556 

2557 

BeO 

2720 

3000  - 

LI 

453 

1620 

li2o 

1843 

2836 

Mg 

922 

1378 

MgO 

3098 

3533 

SI 

1680. 

2630 

S102 

1996 

2500  - 

Zr 

2125 

4776 

Zr02 

2900 

4548  - 

Pe 

1809 

3145 

Pe2°3 

1650 

3687* 

Reference: 

JANAF  (Ref  22) 

♦Decomposes  but  does  not  boll. 
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TABLE  III 


CHEMICAL  ANALYSIS  o.  TRONA  AMORPHOUS  BORON 


Element 

Percer.  ba«e 

(l)* 

(2) 

i* 

Boron 

90.00  - 92.0 

«• 

...... 

Mg 

3.00  - 5.0 

1.7 

- 

3.3 

Fe 

0.20 

0.7 

- 

0.33 

N 

0.10 

- - - 

an 

- - - 

Mn 

0.08 

0.02 

- 

0.04 

Ca 

0.04 

0.07 

- 

0.13 

Si 

0.04 

0.1 

- 

o.3 

Na 

0.04 

- - - 

- 

- - - 

A1 

0.01 

0.007 

- 

0.013 

Pb  . . 

0.01 

0.02 

- 

0.04 

Ni 

0.01 

- - - 

~ 

- - - 

Zn 

0.06 

- 

0.1 

Cu 

Zr 

0.004 

— 

0.006 

0.001 

Absolute  density  2.38  - 2.4 
(gm/cc) 


Average  Particle  Size  1.0  - 1.75  (3)* 

(M) 


References  s 


(1) *  Hand  Book  of  Metal  Powders.  Edited  by  A.  R.  Poster, 

Reinhold  Pub.  Corp.,  N.  Y.,  1966. 

(2) *  Spectrographlc  Analysis,  by  L.  Roush,  Miami  Univ., 

April  22,  1970. 


(3)*  Stokes  Flotation  Analysis,  by  G.  Menninger,  ARL, 
Wright-Patterson  AFB,  Ohio. 
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TABLE  IV 

EXPERIMENTAL  DATA  AT  ^ = 1.0 
(^**)air  £ ik 


22.8 

22.8 


nn/omfa-sec, 

3.20 

3.20 


8.62 

8.62 


0.1 

0.246 


2100 


(cm"1] 

0.055 


1.9 

35.6 

3.20 

13.4  0.3 

2300 

0.44 

0.304 

1.9 

35.6 

3.20 

13.4  0.28 

2400 

0.39 

0.26 

1.9 

35.6 

3.20 

13.4  0.625 

2530 

0.4 

0.268 

2.54 

40.6 

1.80 

27.25  0.9915 

2940 

0.13 

0.055 

2.54 

50.8 

1.80 

34.1  

2950 

0.166 

0.0715 

2.54 

50.8 

1.80 

34.1  

3090 

0.073 

0.0294 

2.54 

5 0.8 

1.8C 

34.1  

2880 

0.153 

0.0655 

2.54 

_50.‘8': 

‘ 1.80 

34.I  

3070 

0.077 

0.0315 

2.54 

33.0 

1.80 

22.15  0.94 

2.54 

56.0 

2.70 

25.0  0.845 

2.54 

56.0 

3.60 

18.73  O.89 

2.54 

6l.0 

2.70 

29.25  0.99 

— 

(/>u 

< 1.5  gm/cm2-sec 

Graphite  Combustor) 

• 

2.54 

15.2 

0.99 

18.6  0.12 

2.54 

25.4 

0.99 

31.0  0.23 

____ 

2.54 

61.0 

0.99 

74.5  0.688 

— 

(/*u  < 1.5  gm/cm2-sec 

Zirconia  Combustor) 

2.54 

30.5 

0.99 

37.2  0.916 

■1  w mm  mm 

2.54 

76.2 

0.90 

94.3  1.0 



2.54 

33.0 

1.34 

29.5  0.49 

2.54 

40.6 

0.90 

54.5  0.715 

2.54 

40.6 

1.25 

39.0  0.372 

2.54 

40.6 

1.44 

34.1  0.424 

— 

(2.54 

cm  Square  Zirconia  Combustor) 

2.54 

20.3 

1.41 

14.4  1.0 
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TABLE  VI 


EXPERIMENTAL  DATA  i.  ? <f>  = 


ID  L 
(cm)  (cm) 

1.9  22.8 
1.9  22.8 
1.9  22.8 
1.9  22.8 
1.9  22.8 

1.9  35.6 

1.9  35.6 

1.9  35.6 

1.9  35.6 

1.9  35.6 

1.9  *35.6. 

1.9  35.6 

1.9  35.6 

2.54*  20.3 
1.9  56.0 

1.9  56.0 

1.9  56.0 

1.9  56.0 

1.9  56.0 

1.9  56.0 

1.9  56.0 

2.54  5Q.8 
2.54  50.8 


2.54 

2.54 

2.54 


50.8 

50.8 

50.8 


(/>*)  air 

3n/cm2-sec] 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 
3.20 
3.20 
3.20 
- -3.20 
3.20 
3.20 
3.20 

1.38 

3.20 

3.20 

1.60 

3.20 

3.20 

3.20 

3.20 

1.80 

1.80 

1.80 

1.80 

1.80 


8.62 

8.62 

8.62 

8.62 

8.62 

13.4 

13.4 

13.4 

13.4 

13.4 

13.4 

13.4 

13.4 

20.8 

21.1 

21.1 

42.2 

21.1 

21.1 

21.1 

21.1 

34.1 

34.1 

34.1 

34.1 

34.1 


0.488 


0.722 

0.77 


0.96 

0.895 

0.973 

0.99 

0.905 

0.89 

0.977 


2250 

2220 

2140 

2250 

2240 

2300 

2320 

2320 

2300 

2250 

2330 

2370 

2330 

2600 

2600 

2500 

2780 

2500 

2570 

2570 

2550 

2630 

2600 

2640 

2490 

2545 


0.25 

0.25 

0.52 

0.58 

0.425 

0.465 

0.265 

0.48 

0.46 

0.38 

0.21 

0.36 

0.36 

0.2 

0.28 

0.165 

0.1 

0.27 

0.29 

0.24 

0.22 

0.24 

0.32 

0.22 

0.22 

0.21 


(cnT*) 

0.151 

0.151 

0.386 

0.455 

0.29 

0.329 

0.162 

0.343 

0.324 

0.25 

0.124 

0.234 

0.234 

0.117 

0.129 

0.095 

0.055 

0.165 

0.18 

0.144 

0.1305 

0.108 

0.152 

0.098 

0.098 

0.093 


* 2.54-om  Square  Zlrconia  Combustor 
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TABLE  VI 
(Continued) 


ID 


2.54 

2.54 

1.9 


2.54 

2.54 

2.54 


2.54 

2.54 

2.54 

2.54 

2.54 

2.54 


L 

(-"“lair 

€ 

V T 

4k 

(cm) 

(pon/cm^-seo)  (ms) 

(°K) 

(cm“l 

30.5 

1.80 

20.5 

0.585  

6l,0 

1.80 

41.0 

0.995  

56.0 

3.20 

21.1 

O.985  

1.5 

gm/cm2-sec 

Graphite 

Combustor) 

15.2 

0.99 

18.6 

0.415  

25.4 

0.99 

31.0 

0.865  

61.0 

0.99 

74.5 

O.936  

: 1.5, 

gm/cm2-sec 

Zirconia 

Combustor) 

30.5' 

1.22 

30.2 

0.945  

30.5 

0.99 

37.2 

1.0  

30.5 

0.99 

37.2 

0.9165  

30.5 

1.18 

31.2 

0.939  

30.5 

0.88 

41.8 

O.993  

30.5 

1.18 

31.2 

0.953  
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TABLE  VII 

TYPICAL  N2/02  AND  Ng/ABGON  DATA* 


L 

(cm) 

* 

V°2 

Ng/Argon 

t 

15.2 

1.0 

4.245 

83.644 

0.12 

15.2 

2.0 

6.205 

84.778 

0.415 

25.4 

1.0 

4.84  5 

85.841 

0.23 

25.4 

2.0 

27.581 

83.57 

0.865 

61. 0 

1.0 

11.988 

83.966 

0.688 

61.0 

2.0  . 

58.681 

84.470 

0.936 

Std.  Air 

\w  ■»  — 

3.728 

83.602 

♦Data  from  2.54-cm  ID  Graphite  Combustor 


147 





v'  ygggg 


-V~\% 


»«w  T’ 1 '^cv^^^^^yy^yyyaTOr*^^ 


| 


{ 


> 

I 

I 

! 

; 


DS/ME/74-2 

VITA 

Lowell  W.  Ormand  was  born  on  2 April  1932  near  Butler 
in  Custer  County,  Oklahoma.  He  graduated  from  Amphitheater 
High  School,  Tucson,  Arizona,  in  May,  1950.  In  May,  1955, 
he  received  a Bachelor  of  Mechanical  Engineering  from  the 
University  of  Arizona  and  received  a commission  as  a second 
lieutenant  in  the  USAF  Reserve.  He  was  employed  by  the 
Shell  Chemical  Company  as  a plant  maintenance  engineer  until 
called  to  active  duty  in  July,  1955.  His  Air  Force  assign- 
ments before .coming  to  the  Air  Force  Institute  of  Technology 
were  as  a combat  crew  member  in  the  303rd  Bomb  Wing,  SAC; 
student  with  the  AFIT  Civilian  Institutions  Program  at  the 
University  of  Arizona  where  he  received  a Master  of  Mechanical 
Engineering  in  May,  1965;  project  engineer  in  the  6511th 
Parachute  Test  Group  and  with  the  7th  Air  Force  in  Vietnam 
as  a forward  air  controller.  He  was  assigned  to  the  Air 
Force  Aero  Propulsion  Laboratory  where  he  accomplished  the 
present  investigation  and  is  currently  assigned  to  the  Air 
Force  Flight  Test  Center.  He  attended  pilot  training  in 
1955*  Squadron  Officers  School  in  1959  and  accepted  a 
commission  in  the  Regular  Air  Force  in  1962. 

Colonel  Ormand  is  married  to  the  former  Barbara  Welch 
of  Tucson,  Arizona  and  they  have  four  children,  David, 
Katherine,  Michael  and  Rebecca. 

Permanent  Address:  2227  E.  Drachman 

Tucson,  Arizona  85719 

This  dissertation  was  typed  by  Barbara  Ormand; 


\ 


148 


